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1 Introduction 

Over the past few years there has been growing interest in using gauge/gravity duality to 
model strongly interacting systems relevant for condensed matter physics, see the reviews 
[1-4]. In particular, we have learnt that AdS black holes generically develop hair at low 
temperature, with the hair capturing the phenomenology of superfluid phases in the dual 
field theory [5, 6]. 

It is remarkable that rather simple and generic gravity models can capture many features of 
the phase structure of strongly interacting systems but to construct more sophisticated and 
more realistic models one clearly needs to include additional ingredients. Every essential 
ingredient or feature of the strongly interacting system must have a counterpart in the 
bulk holographic model. There has been considerable interest recently in developing bulk 
descriptions for key missing ingredients such as broken translation invariance (necessary to 
realise finite DC conductivity in the ordinary state) [7, 8]. 

Moreover one needs to take into account that there can be a rich variety of phases even 
within rather simple models. Isotropic superfluid phases are not necessarily the preferred 
phase at low temperature: spatially modulated phases arise in rather generic models, see 
for example [9-13]. One may also find other phases, including striped phases [14], helical 
phases [15], and simple models can exhibit competitions between such phases [16]. 

The focus of this paper is on realising an important missing phase, d-wave superfluidity, 
within a top down holographic model. The importance of modelling d-wave superfluidity 
is self-evident: many unconventional superconductors admit either d-wave or mixed sym- 
metry. A natural candidate for modelling the d-wave condensate is to use a charged spin 
two fleld in the bulk, instead of a charged scalar field. This approach was taken in [17-21] 
and also in [22, 23]. The problem is that massive spin two equations have issues with con- 
sistency, causality and hyperbolicity, and they do not follow from consistent truncations of 
reductions of top down models. 

The approach in [17] and subsequent related works was to write down minimal equations 
for the spin two field (without looking in detail at the constraint equations required to 



obtain the correct number of propagating degrees of freedom) while [22, 23] looked in more 
detail at the effective action for the spin two field and how the constraint equations could 
be satisfied. In practice, evaluated on a static ansatz, the effective equations of motion 
are the same in both approaches. The coupled equations reduce to those for the charged 
scalar field/gauge field system used in [6] and therefore an analogous condensed phase 
exists at low temperature. Evaluating the onshell action, one can show that it reduces 
to the same action as for the Maxwell-scalar case, and therefore the condensed phase is 
thermodynamically favoured. 

The starting point of this work is the observation that massive spin two modes are generic 
in compactifications from ten and eleven dimensions which give rise to lower-dimensional 
Anti-de Sitter backgrounds. Moreover such modes are generically charged under the C/(l) 
symmetry corresponding to the R symmetry of the dual conformal field theory and thus the 
corresponding spin two operators in the dual CFT are natural candidates for d-wave order 
parameters. (An analogous observation was made for circle reductions of geons in [24], 
although the charged spin-2 condensates obtained in this case were not thermodynamically 
preferred.) 

While charged spin two fields are generic, they cannot be retained as part of a consistent 
truncation and therefore one has to derive their field equations directly from the higher 
dimensional equations of motion. The linear field equations suffice to obtain the spectrum 
but, just as for the coupled system of a charged scalar and gauge field, non-linear inter- 
actions are essential in [17, 23] for the condensate phase to exist. The problem is that 
as soon as one allows interactions between Kaluza-Klein modes Pandora's box has been 
opened and one does not expect to be able to truncate to a finite set of equations: the spin 
two field will source an infinite number of other Kaluza-Klein fields. 

Clearly the system of equations can only be tractable if one can suppress the backreaction of 
the spin two field on the other Kaluza-Klein fields, including the metric. There are (at least) 
two distinct possibilities for suppressing the backreaction. The first is to work in a limit 
in which the charge q of the spin two field is large, following the arguments given in [4, 6] 
for a charged scalar field, see also the discussions around (2.48) and (2.49); then the idea 
is that the backreaction would be suppressed by 1/q factors. The second possibility is to 
work close to the critical temperature Tc when the spin two field magnitude is proportional 
to e = (1 — T/Tc)^ <C 1 (for a specific exponent A); the backreaction of the spin two is then 
suppressed by factors of e^ or smaller. In this paper we will consider both possibilities but 
the analysis is simpler and more robust in the latter case. 

Large charge q certainly suppresses the backreaction, at least in part. However, in a 
top down model one needs to take into account that the mass and the charge are not 
independent: there is a maximal charge for any given mass. In particular this implies that 
the critical temperature is not reached for spherical or simple Sasaki-Eintein reductions 
within the approximation that the backreaction on the metric is ignored. (The probe 
approximation may however suffice for other reductions giving AdS vacua, such as warped 
compactifications.) The issue is that the backreaction of the gauge field on the metric is 



controlled by /x/T where // is the chemical potential and is only negligible when ^/T <^ 1. 
The mass and charge parameters in the compactifications are however such that Tc/fj, < 1. 
This implies that one must take into account the backreaction of the gauge field, which in 
the ordinary phase gives rise to an R charged black hole. 

Suppose one considers a spin two field in the (fixed) background of an R charged black hole. 
Looking at the linear equation for the spin two field in such a background, one finds that 
a condensed phase indeed generically exists below a critical temperature. One does not 
need non-linear equations for the condensed phase to exist. This might sound surprising 
given the discussion above but in fact this is exactly what happens close to the critical 
temperature in the s-wave case: the non-linear coupled equations of the gauge field and 
charged field are solved perturbatively in the parameter e given above, and the leading 
order equation for the charged field becomes a decoupled linear equation in this limit, see 
[25-27] and appendix C. 

The linear equation for the spin two field is however only a valid approximation if the 
backreaction of the spin two onto the metric, gauge field and other Kaluza-Klein fields is 
suppressed. The backreaction onto the metric and the gauge field is suppressed either for 
large q (with the amplitude of the spin two field scaling as 1/q) or close to the transition 
temperature. The backreaction on the other Kaluza-Klein fields is however more subtle 
because it involves looking in detail at the interactions in the Kaluza-Klein reduction. Here 
we will argue in detail that the other Kaluza-Klein fields can indeed be neglected close to 
the critical temperature, postponing a detailed analysis of the large q case for future work. 

Close to the critical temperature the description of the spin two condensed phase is sim- 
ple: a linear perturbation around the charged black hole background. This description can 
straightforwardly be uplifted to ten or eleven dimensional supergravity; the spin two field 
then corresponds to a specific graviton perturbation in the background of decoupled spin- 
ning D3-branes or M2-branes. This upstairs picture would be the natural starting point 
for the construction of a fully non-linear condensed phase solution. 

Although the focus of this work is on d-wave phases, Kaluza-Klein modes may be re- 
sponsible for a rich variety of phases, many of which may not be visible within consistent 
truncations. The d-wave phases we find would have critical temperatures which are typi- 
cally lower than those of s-wave phases associated with relevant scalar charged operators. 
In a compactification in which no such relevant operators are present, the d-wave phase 
may be the first superfluid phase reached as one lowers the temperature. In typical com- 
pactifications, the s-wave phase would however be reached first (unless one finds a way to 
lift the scalar modes). One would expect that the interactions between distinct Kaluza- 
Klein modes would give rise to a competition between different phases, which would be 
very interesting but perhaps complicated to analyse. 

In the bottom up model, the onshell action for the charged spin two field system is com- 
puted, showing that the condensed phase is indeed favoured below the critical temperature. 
Moreover, in the probe approximation, spin two and scalar fields of the same mass and 
charge would have exactly the same free energy in the condensed phase, neglecting any in- 



teractions between them. Taking into account either interactions or the backreaction onto 
the metric would presumably break the degeneracy and determine what is the favoured 
phase. 

Working near the critical temperature in the top down model, one cannot determine 
whether the condensed phase is indeed thermodynamically favoured without computing 
the backreaction, both on the metric and on the other fields. The most straightforward 
way to compute this backreaction would be to work directly with the higher dimensional 
description but we postpone such an analysis for future work. 

This paper is organised as follows. In the next section we consider the effective description 
of d-wave condensates in terms of coupled spin two and gauge field equations in black hole 
backgrounds. In section 3 we discuss how such equations may arise from Kaluza-Klein 
reductions and in section 4 we look in detail at reductions of type IIB on Sasaki-Einstein 
manifolds. In section 5 we present the description of the condensed phases in terms of 
higher dimensional solutions and in section 6 we conclude. Various technical issues are 
contained in the appendices: appendix A relates to corrections to the spin two equations 
in Kaluza-Klein reductions while appendices B, C and D explore numerical and analytical 
solutions of the condensate equations. 

2 Effective description of (i-wave condensates 

We will consider holographic superconductors in {D + 1) spacetime dimensions, using an 
action involving the metric g^y together with a spin two field 0^jjy and a gauge field F^y. 

There has been considerable literature discussing actions for spin two fields coupled to 
electromagnetic fields. Working from a bottom up perspective, [28-31] discuss the con- 
ditions under which such fields can have consistent and causal equations. Inconsistencies 
and acausalities are analysed in [32-34] . Here we will use the equations of motion following 
from the action given in [22, 23]^ 

S = ^ f d^^'xV^ (r-2A- \f^uF^^A (2.1) 
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2^2 



+ ^_ / d^ + l^;, 



\Dpct>^y\^ + 2\D^r''\^ + \D^^\^ - [D^<lf^''Dy(t) + c.c] 



1 

D + 1 



-m" (|</>^,|2 - |,^|2) + 2R^ypx4>*^P4>''^ - j—rrR\ct>? - 2igqF^y4>*^^4>l 



Here (f) = (j)^ and (f)p = D^^cpf^p, with the gauge covariant derivative being defined as 
D^ = V^ — iqAp. The parameter g is the gyromagnetic ratio which will be discussed 
further below. The Chern-Simons terms involving the gauge field which are present in odd 



^Note that there is a typo in the Lagrangian given in equation (3) of [23]; the term involving Rii^cj>*'^'''4>'\ 
should be absent, following the discussions in appendix A of [23]. The equations of motion in equation (4) 
of [23] indeed drop this term. 



dimensions do not play a role in what follows, as we consider only static electrically charged 
solutions. 

This action admits Einstein solutions in which the spin two and gauge field vanish, and 

the metric satisfies 

2A 

^^ ~ (D- 1) ^^"" '^^■^'* 

We will fix A = —D[D — l)/2 so that the Einstein metric has unit radius. For much of 

this paper we will treat the spin two and gauge field in the probe approximation, in which 

their equations of motion are solved in a fixed background Einstein metric. The equations 

of motion then describe a massive charged spin two particle of mass m and charge q in the 

curved Einstein background: 

(□ - rr?)(t>^y = 2L><^0^> - D^^D^yc/) + g^^ [(D - m^)(t) - DP(t)p] (2.3) 

n 

+2R^pXu4>''^ + 9,.u-rjj-p^^ + igq {F^p€ + ^^pO ; 
Z^^F^^ = iqct>% (^D'^cP^f - 2gD^cP^P 

-iq(2{l - g)cl>;rP - D,cl>*r' + (I?>* - <Pn 0j + h-c 

(Anti)-symmetrization (without subtracting the trace) is denoted by angled parentheses 
and is defined with unit weight. These equations of motion are supplemented by constraints, 
which when F^i, = take the simple form 

D^cP^, = = 0. (2.4) 

In this case the equation of motion reduces to 

(D - m^)(l)^, + 2R^xvp4>^' = 0, (2.5) 

subject to these constraints. These constraints arise as follows: the components of the 
equations of motion in the time direction do not involve second order time derivatives and 
are therefore constraints. The tensor cp^n has more degrees of freedom than a spin two 
particle, and additional constraint equations are needed to eliminate the extra degrees of 
freedom. These further constraints arise from taking the divergence of the equations of 
motion twice. 

In the case where F^y ^ the constraints are considerably more complicated. Only the 
case oi g = 1/2 is discussed in [23] but we will need the equations for general g which were 
given in [33] . The first divergence of the field equation gives 

w?{(t>^ - D^(t>) = -iq Qf^,^ - (1 - g)F^(t>u^. + gOpF^^cf"' (2.6) 

(1 + g)Fp,r - \Fp,D^ct> + (2 - g)F''PDpcPp, 



where F^ = D^F^^. The trace of the field equation gives 

D + l 



Dn?--^^—\R](t)-{D-l){n(t)-Dp(P^') = {). (2.7) 



The second divergence of the field equation combined with the trace equation gives 
Dm'^ R 



m 



{D-l) {D-l) 



= iq [2(1 - g){D^F^P){D,(l,p,,) - IgF^^cP^ + F^D^cP] 

+ {g - 2)q-'F'^PF,pr, + -/F^,F^^'^ct> (2.8) 

+iq{l - 2g)F^''D^cl), + iq{l - 2g)DpF,cl)^'\ 

From these equations it is apparent that when ^^ = the constraints reduce to (j) = (f)^ = 0. 

When Ap ^ 0, the problem is that these constraints are no longer algebraic constraints on 
{(j),<j)f^); they involve first derivatives of (p^i, multiplied by first derivatives of F^^i, and for 
generic values of F^i, this leads to the equations being either non-hyperbolic or non-causal 
[33]. Prom the form of (2.8) it is apparent that for generic values of g there is a term 
multiplying the field strength with double time derivatives of the spin two field. Such a 
term implies a breakdown in the number of degrees of freedom, as a constraint of the free 
model becomes a propagating degree of freedom. This argument has been given in the 
literature, see [33], [23], to fix ^r = 1/2, in which case the action describes (-D + 1){D — l)/2 
propagating spin two degrees of freedom. 

As we see below, the equations and constraints can be consistently solved for specific static 
field configurations, such as those discussed below, for any value of g. If one then analyzes 
perturbations around such a solution, generically one will have a number of degrees of 
freedom which does not match that of a spin two field. 

Even for g = 1/2 in a fixed Einstein background, the equations of motion for general values 
of the gauge field are still non-hyperbolic or lead to non-causal propagation. To suppress 
these effects one can work in a limit in which 



(7 F-- n\F~~ - 



(2. 



are small, where fi, denotes a tangent space index. It was suggested in [22, 23] that the 
problems with causality and hyperbolicity could be corrected by adding additional higher 
order terms to the action. 

Most of the discussions of spin two and gauge field actions in the literature follow from 
bottom up models. Here we are working in the context of Kaluza-Klein reductions of string 
theory/M theory backgrounds compactified on Sasaki-Einstein spaces. Kaluza-Klein gravi- 
tons in such reductions correspond to massive spin two fields from the lower dimensional 
perspective. In the following sections we will show that there is indeed a systematic way 
of expanding the spin two/gauge field equations when they are obtained from a Kaluza- 
Klein reduction. In anticipation of this discussion, it is useful to note the following point. 
Suppose that one treats both the gauge field and the spin two field perturbatively, i.e. one 
lets 

<^/.. = S^i^ + <5'</'?j + • • • ; A^ = M« + J^^if) + • • • , (2.10) 

with 5 small. Working perturbatively in 5, the constraint (2.8) remains algebraic. To 
leading order we obtain that the trace i;^'^' = and to next order we obtain a constraint 



on (/)(^) in terms of the (known) fields (fy^l and A^ , i.e. we obtain from (2.8) 

- ""' OD^ (^' + (^ + ^^) '^^'^ = ^^^(^ " 9){D>^F'^^^'''D,^^pi)- (2-11) 

Thus if the equations are evaluated perturbatively the restriction to g = 1/2 is not neces- 
sary. From a top down perspective the parameter g is uniquely determined by reducing the 
higher dimensional equations of motion over a compact space. Constraints on such cou- 
plings may be obtained by imposing causality in the dual CFT, see [35]. The top down will 
automatically respect such constraints but furthermore determine such couplings uniquely. 

2.1 The d-w^acve condensate 

Let us now turn to superfluid solutions of these equations of motion. The Einstein condition 
on the metric restricts us to considering a probe limit in which the matter fields do not 
backreact on the metric. The relevant background metric is therefore that of a neutral 
black brane 

ds^ = 1 (-fiz)dt'' + dx^ + ^] , (2.12) 

where 

/(.) = l-(i)". (2.13) 

The black hole horizon is located at z = Zh whilst the conformal boundary of the spacetime 
is at z = and the temperature of the black hole is 

T= . (2.14) 

A-KZh 

We now consider solutions to the coupled equations of motion for the spin two field and 
gauge field, assuming an ansatz in which cp^^i, and A^ depend only on the radial coordinate 
z, with only spatial components of the spin two field switched on. Furthermore, we consider 
the case in which only (p^y 7^ is non-zero, so the ansatz is 

A = V{z)dt- 4>xy = J^Kz), (2.15) 

with the functions {V{z),h{z)) being real. Note that this is clearly not the most general 
ansatz consistent with the symmetries: one could obtain related solutions by imposing a 
complex ansatz for the spin two field and by breaking the rotational symmetry in a different 
way. Put differently, there is an arbitrariness in the choice of the phase characterising the 
[7(1 ) symmetry breaking and in the choice of the plane in which the rotational symmetry 
is broken. With this ansatz (f) = (f)^ = F^p(p^ = and the equations of motion simplify 
considerably - the only terms which contribute are the following 

(D - m^),^^, = 2R^xp.r^; (2.16) 

D^F'- = igC«(D>-^)+h.c., 



where □ = D^D^. Note that the gyromagnetic parameter g drops out of the gauge field 
equation, because of the symmetries of the background and of the ansatz. One also needs 
to check that the constraints are satisfied. The conditions cj) = cj)^ = F^ycf)^ = are already 
sufficient for (2.8) to be satisfied but the remaining constraints (2.6) and (2.7) require the 
stronger conditions of (2.15). 

Before moving on to discuss these solutions of the equations, we should note that the 
coupled system could of course admit many more types of solutions and hence different 
phases. In what follows we will concentrate on the simplest d-wave superfluid phase, and 
we will postpone discussions of other phases and indeed the related question of the stability 
of this superfluid phase for later work. 

Equations of motion 

Using the ansatz (2.15) the coupled equations of motion become 

V" + ^^J^V' - t^h^V = 0; (2.17) 

where a prime denotes a z derivative. These equations of motion admit a normal phase 
solution in which the spin two field is zero 

dt- (P^y = 0, (2.18) 

with /i being the chemical potential. 

In [22, 23] it was noted that these equations of motion are identical to those of the Abelian 
Higgs model in AdSo+i discussed in [6], implying that there also exists a condensed solution 
below a critical temperature Tc. In the condensed solution the spin two field is non-zero, 
and its asymptotic expansion is 




K 



2 



h{z) = z^-^ {hs + 0{z)) + z^ (^^^—^(O,^) + 0(z)J , (2.19) 

where A is the conformal dimension of the dual spin two operator and the source term h^ 
vanishes. The conformal dimension of the spin two field is related to the mass parameter 
rm? as w? = A(A — D). We will derive this relation between the normalisable mode and 
the expectation value of the dual operator {Oxy) in section 2.5. Moreover, in the condensed 
phase regularity conditions are imposed in which h{z) is finite everywhere and the gauge 
field At vanishes on the black hole horizon. Numerical solutions for the condensed solution 
in D = 3 were presented in [22, 23]. Note that the non-linear interactions in the equations 
(2.17) are necessary for a condensate solution to exist. 

For later use, let us notice that the ansatz is such that only two terms non-linear in 
{(pxy^^t) appear in the equations of motion in (2.16). These non-linear terms follow from 
the covariant derivatives implicit in the action (2.1). 



We now introduce a dimensionless radial coordinate C = z/ Zh such that <^ = 1 at the 
horizon and C — ^ at the boundary. Then let the gauge field and spin two field be rescaled 
as 

V{C) = qzhV{0- h{C) = qh{C). (2.20) 

The coupled equations of motion then become 

n • 1 

(2.21) 




with f{() = (1 — (^) and d = d(^a. The rescaled equations depend on only one param- 
eter, the mass of the spin two field, and are independent of both the charge q and the 
gyromagnetic parameter g. At first sight this implies that the existence of a non-trivial 
condensate solution depends only on the mass; however, the rescaled equations need to be 
solved subject to the boundary conditions 

V{C = 0)^f^^ mo; V{C = 1) = 0, (2.22) 

47ri 

with the spin two field h having no source term as (" — ^ and being finite at the horizon ^ = 
1. This implies that the existence of a condensate solution does indeed depend implicitly 
on the temperature, charge and chemical potential in the combination mp. For 

nio > acin (2.23) 

with Oc a (numerically) determined constant the condensate solution will exist and be 
preferred. The parameter Qc corresponds to a critical temperature 

- ""^^ (2.24) 



47ra^A 



with A the dimension. Note that etc ~ 3.68(2.41) for D = 3(4) for large A ~ m. For small 
values of m in Z? = 3 see Table 1 . For D = A the parameter Uc also increases monotonically 
with dimension and is bounded as Oc > 1.4, with the lowest value obtained at the lowest 
dimension, A = 2. 

In the next subsections we consider how this critical temperature may be estimated for small 
and large values of m. The reason for analysing the temperature is that in top down models 
a diverse range of the parameters (g, m) may be realised but they are not independent of 
each other. It is therefore useful to study in detail how the critical temperature varies as 
these parameters are adjusted, without solving the equations numerically, case by case. 
Note that numerical solutions become increasingly difficult as one increases m. 

Estimation of critical temperature 

In the case in which the dimension of the spin two operator is large one can estimate the 
critical temperature following the arguments given in [22, 23]. The idea is to consider a 



A+K) 


2(-2) 


3(0) 


3.5(7/4) 


4(4) 


5.25(28) 


oo 


ac 


2.07 


2.52 


2.69 


2.806 


3.28 


3.68 



Table 1. a^ as a function of A+ from numerical analysis [D — i). We substitute the dual operator 
dimension A+ for m in (2.23); for relevant operators we consider the quantisation leading to a 
higher dimension, hence the notation A-|_. 



massive charged point particle in the black brane background in the ordinary phase, in 
which there is a background gauge field. For mass m and charge q the point particle action 
is 

S = -m I dr-q I A. (2.25) 



Then the effective potential for a static particle is given by 



Kff(C) 



m 



D~2\ 



(i-n^^ + gMi-C""^) 



(2.26) 



with ^ = zjzh. The extremum of the potential is at V^g(Co) = 0; this will be a global 
minimum only if q\x is large enough such that 14fr(Co) ^ 0, since the potential is zero at 
the horizon C, = \. In D = 3 both conditions are satisfied when [22, 23] 



\qpL.Zh\ ~ 3.68m, 



(2.27) 



which agrees well with numerical results (but slightly overestimates the critical tempera- 
ture). In D = 4 the resulting condition is rather similar: 



n^zh\ 



which corresponds to a critical temperature 



2.41m, 



tr - 7.57A, 



(2.28) 



(2.29) 



where we use the fact that A ~ ttt, for large dimension. The corresponding value of Co is 
found to be 0.64, which is outside the horizon. Condensation thus occurs at temperatures 
which are small compared to q\i\ there is no condensation when the chemical potential is 
switched off. 

The actual value of the condensate can be obtained by numerical solution of the field equa- 
tions, with the charge density p and the chemical potential \x read off from the asymptotic 
values of the gauge field potential 



Y[z) = ii- pz 



D-2 



(2.30) 



For example, in the case of Z? = 4 with rr? = 5, A = 5, 5 = 1, figure 1 shows numerical 
condensates for which K^{Oxy) ~ (p/p) '^^~'^' with a coefficient of order one below Tc. 
(Note however that in this case Tc/p is not small, and therefore the probe approximation 
is not justified. We will return to this issue below.) 
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Figure 1. Condensate: £) = 4,771^ = 5, A = 5, g = l,c = ^"^/'^(p/Ai)^/^^"^). 



2.2 Analytic approximation 

Consider the dimensionless equations (2.21). In the near boundary region, (^ — )• 0, and in 
the near horizon region, around C = Ij the solutions take the form 



^(0 



MO 



Vh = vo- vi 



7./-D-2 



if C^O 
(l-C) + i(-3 + I?-|)(l-C)'l ifC^l 



/lB = /l+C^ + /l_C 



d-A 



hu = h, 



1 + 1^(1-0 + (^^-+--"^) (i-C)^ 



if C^O 

ifC^i 



(2.31) 



The near boundary data, /z, p and /i+, are proportional to the chemical potential, the 
charge density, and the one point function of the spin two operator respectively: 

A 



,47rr 



P = P' 



qD q\Dj K^V^q \ D 



P = P- 



(2.32) 



The coefficient /i_ corresponds to the source of the operator and thus we set /i_ = 0. 
Imposing regularity at the horizon, we choose uq = for V and have dropped the other 
independent singular solution of h. 

Now let us match the solutions at an intermediate value ^m, requiring continuity of the 
functions {V , h) and their first derivatives. This gives four conditions 



VBiQm) = VHiCm) , ^b(U) = V H{Cm) , 



(2.33) 
(2.34) 



with six unknowns: p, p, vi, h+, ho, Cm- We will express four variables p, vi, h+, ho in terms 
of p, Cm- First, from 

hBiCm) hniCm) 



hB\Qm) il-HyCn 



we have vi as a function of Cm only 



o o 4 2DA(m'^ + D(2-A)) 2D(m'^ + DA) 



A(l - Cm) + 2C« 



l-Cv 



(2.35) 



(2.36) 
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The second equation of (2.33) yields 



^"1(1 -Cm)/D 

with which, the first equation of (2.33) gives 

- Via -Cm) ^ pC-HD{l-Cm) + ^Cm-2) 

fl - \ . (2.38J 

Before discussing h+, let us fix the sign of vi and /iq- The expression of /iq (2.37) gives 
some constraints. First of all, [4 — D + (D — 3)Crn] is mostly negative for D > 5. However, 
if one chooses r]m such that [A — D + {D — 3)Cm] < 0, we cannot see the condensed phase, 
whatever the sign of vi is. So, this physical consideration restricts ourselves to D = 3, 4 and 
[4 — D + (D — 3)(^m] > 0. This restriction may however be due to the limitations of our crude 
approximation. The second constraint is that only vi > is allowed for h^ > 0. Thirdly, 
the sign of /iq is arbitrary and does not matter by symmetry so we choose a positive sign. 

Finally, /i+(~ {Oxy)) is solved by (2.34) with (2.36) and (2.37). It is written in terms of 
physical variables as 



2A-_D+1 / / rr\ D-l 



^O.y) _ ,,^^ , ._ „^./rT^— \/^ (T 



^^ -A{D,m,A,K,Cm,q)^[l^j \l^-[Trl ' ^^-^^^ 



where 



2A-D /47r\^ (2Z) + m2(l-U))C^ D{4 - D + jD - 3)Cm) . . 

K^V^q\Dj DA-D{A-2)Cm V ^ " Cm ' 

"^^ = i^ Ui4-D + iD-3)u) ' ' ^'-''^ 



I ^ , 2D(-m^ + D(A-2))A 2DA^ , , , 

vi = J 2Dm^ + m'^ ^-— -—^ t^^ + :- ■ (2.42) 

V A(l-Cm) + 2Cm I -Cm ^ ' 

A real solution for (Oxy) exists only when T < Tc and the critical exponent is 1/2 for all 
parameter set {Oxy ~ \/l — T/Tc). 

With a matching point ^^ = 0.5, the formula is simplified as 



■^~ K^V^qKo) A + 2 ' ^ '^> 

llGD^A + m4(2 + A) + 2Dm?{<o + 5A) 

^^ ~ V 2 + A ■ ^ ■ -^ 
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For example, with the parameters {D = 3,m? = —2, A = 2,-^/2^ = l,q = 1), (2.39) is 
reduced to 

/m \ 1 R„2 O I /- I 5 rp I / rp \ 2 

(2.46) 

(2.47) 

which agrees with (A. 30) and (A. 31) in [36]. (Note that such a parameter choice is vahd for 
scalar fields but is invalid for spin 2 fields, as m? > for unitarity.) One has to treat this 
approximation with some caution since the matching point was arbitrary. However, the 
simple analytic calculation is in good agreement with the numerical results, particularly 
close to the transition temperature. Further discussion of the structure of the coupled 
equations is given in the appendices and will be used in the later discussions. 

2.3 Validity of approximations 

Now let us turn to the conditions under which the use of the probe equations of motion 
can be valid. There are two main issues. First of all, the backreaction of the matter fields 
on the metric has been neglected. Secondly, the action is not causal, although to suppress 
acausal effects one can require that all components of the gauge field strength in a tangent 
frame are small relative to rr? jq. 

Looking at the first issue, notice that in the ordinary state (in which the spin two field 
is switched off) this condition is actually unnecessary as Einstein gravity with negative 
cosmological constant and gauge field is a well posed theory, which can moreover be ob- 
tained as a consistent truncation of supergravity reduced on a Sasaki-Einstein manifold. If 
nonetheless we consider the limit in which the backreaction of the gauge field strength on 
the metric can be neglected in the ordinary state, then necessarily 

T~» ix. (2.48) 

Noting that the critical temperature T^ is smaller than q\i the neglect of the backreaction 
is only consistent with reaching the condensate phase if fi <^Tc < qfi, which requires large 
q. This is the standard limit in which the probe approximation applies, following [4, 6]. 
More precisely, however, looking at the condition for the critical temperature (2.24) one 
sees that the condition required is 

q » — ^. (2.49) 

In a bottom up model in which q and m are independent, one can always choose q such 
that this condition holds. In a top down model, or indeed taking into account the bounds 
on R charge for operators of a given dimension, q cannot be arbitrarily large for a given 
m. In the following section we will consider the case where this large charge condition is 
dropped. 



13 



Let us turn now to the second issue of causality and hyperbolicity violations. It was pointed 
out in [22, 23] that these could be suppressed by imposing (2.9), which in turn requires 
that the dimension of the spin two operator is such that A ^ 1. One should note however 
that such a condition is not needed for the static condensate solution to exist; it would be 
needed in analysing the dynamics of fluctuations around the condensate. In other words, for 
the effective action to be valid for discussing the stability of the condensate, perturbations 
around the condensate should behave causally. If they do not, then the action is inadequate 
for describing the stability of the condensate. However, in the later sections of the paper 
we will argue that in top down models one will in any case need to take into account other 
fields when discussing the stability, and the stability is best analysed directly from the 
higher dimensional equations of motion. 

2.4 Near critical temperature 

We now consider the case where the parameter choices {q, m) are such that the critical 
temperature obtained is not large compared to the chemical potential. In this case the 
probe approximation is insufficient to analyse the existence of a condensed phase. However, 
one can analyse such cases as follows. In the ordinary phase the backreacted (Reissner- 
Nordstrom) metric can be written as 



ds^ = \ (-F{z)de + dx^ + -^ ) ; (2.50) 



,,2 2 D 

F{z) = f{z) + ^'^' 



where f{z) was given in (2.13) and 




9 2(D-1) , , 



The temperature is now given by 
with 



D ( iD-2Wzl \_ D ( (Z^-2) \ 

V = ^, (2.53) 

T 

where Zh still denotes the horizon location and the parameter < f] < D/{D — 2). One 

can also write 

T D / (D-2) \ , , 

'l-^--^ry . (2.54) 



fi 47r7y^ \ D 

The extremal limit of the black hole is obtained when r] — > D /{D — 2) with the uncharged 
limit being ry — )• 0. 
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We now consider the static ansatz for the spin two field, (pxy = l/(2-2^)/i(-?), in this back- 
ground. Close to a (continuous) phase transition the condensate field h{z) would scale 
as e = (1 — T /Tc)^ <^ 1 where A determines the critical exponent, see the discussions in 
appendix C. The exponent obtained from the non-linear equations is A = 1/2. In this limit 
the backreaction of the spin two field on both the metric and on the gauge field would be of 
order e^ and thus negligible. So to analyse the existence of a condensate phase we need to 
look at the linearized equation of motion for the spin two field in the charged black brane 
background. Using the spin two equation given in (2.3) for a transverse and traceless field 
we obtain 

(n - m^) ^^, = 2R^pxA'^. (2.55) 

The resulting equation is very similar to that obtained in (2.17), 

Rewriting in terms of the dimensionless variable C, one obtains the linear equation 

^ ^-^^■- ^^^ "^'\ = 0, (2.57) 





with F(C) = (1 - C^ + r7C^(C^-2 - 1))) and 

V = qfiZh (1 - C"""') • (2.58) 

The linearised equation needs to be solved subject to the condition that the field is regular 
at the horizon and has no source behaviour at infinity. Prom the rescaling it is manifest that 
the condensation condition can involve only the parameters {m,mo = q^Zh,Zh) together 
with the dimensionless number rj which characterises the non-extremality of the black hole. 

The critical temperature Tc can be parameterised as 

T. = -^-^(l-(Z?-2)^). (2.59) 

where the parameter ac{rj,ra) now depends on 77 as well as on m. For fixed chemical 
potential one can combine this equation with (2.54) to obtain a bound on the parameter r/ 
for a condensate to form: 

Tl> 2 2- (2-60 

Repeating the estimates of the critical temperature and the analytic approximations, one 
finds that adrj) is actually only weakly dependent on rj. For example, if one uses the 
effective potential technique for high masses, the relevant potential is now 

Kff oc i (1 - C^ + r?C^(C^-2 - 1))^/' - ^ (1 - C^-2) (2.61) 

(^ ' \m\ 

with OLcij]) being the value of \molvn\ for which the potential admits a minimum at a 
negative value. The parameter Oc decreases with increasing r] but remains of the same 
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order. In Z) = 3 at large m we find adO) ~ 3.68 and ac(1.5) ~ 3.35 while ac(3) ~ 2.46 
(with r/ — 7- 3 being the extremal limit). In Z) = 4 at large m we find ac(0) ~ 2.41 and 
ac(l) ~ 2.24 with ac{2) ~ 1.74 (where r/ — t- 2 is the extremal limit). 

Returning to (2.60), and recalling that r] < D/{D — 2), we see that a condensate only forms 
at finite temperature if the mass and charge parameters of the condensate satisfy 

m2 2D{D - 1) 

^ - (D - 2)2^2 (2-62) 

For Z? = 3 a conservative estimate (at large m,) would be rr? jq^ < 4 and in Z? = 4 at large 
m we find that m? /cf' < 2. 

For smaller values of the mass we can use the methods discussed in appendix C.2 to 
estimate Oc- As one would anticipate from the previous results, Oc decreases with decreasing 
operator dimension and decreasing temperature, i.e. closer to the extremal limit. Looking 
at the lowest operator dimension in Z) = 2, A = 2, one can estimate that ac(0) ~ 1.4 while 
ac(1.75) ~ 1.1. (For larger values of r] the approximation becomes increasingly sensitive to 
the details of the stability criteria.) Therefore in the low dimension limit in Z) = 4 we see 

that 

\m \ 

^ < 6, 2.63 

where the magnitude is used since for relevant dimension operators m? < 0. A similar 
bound was discussed in [4]. We will use these facts later when discussing the modes which 
arise in Kaluza-Klein reductions. 

2.5 Thermodynamics 

The existence of a condensate solution does not guarantee that the condensate phase is 
thermodynamically favoured over the ordinary phase. In this section we will analyse the 
thermodynamics and show by computing the onshell action that the condensate phase is 
indeed preferred. 

Naively one might think that one cannot compute the onshell action to the required order 
without analysing the backreaction of the metric. In the probe approximation, we work 
with a fixed background which solves the field equations and solve the gauge field/charged 
field equations in this background, for small chemical potential fi <^T and fixed {m, q). If 
one looks at the contribution to the onshell action from the gauge field, it would appear to 
be of the same order (/i^/T^) as the backreaction of the gauge field and charged field on 
the metric, which suggests that one needs to know the backreaction to compute the onshell 
action. 

This naive reasoning is not correct, however. To understand this point let us consider 
a more general situation where one has an action S[g, $ ] involving a metric g^i, and 
a collection of matter fields $ (which are not necessarily scalar fields). Let g° and 
($°) be a solution of the equations of motion following from the action and consider 
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perturbations around this solution such that 



M 



-^o\M 



($«)JW + A$ 



M 



The onshell action for the perturbed solution is then 

5'onshell(5) ^ ) = 'S'onshcU b° , {^°)"^] + (^S'onsholl 

with the change in the action being computed using 

SS = I d'^+^x 






6g' 



■ flU 



§^M 



(2.64) 



(2.65) 



(2.66) 



with C being the Lagrangian. By definition the action is stationary on any onshell solu- 
tion. In particular, this implies that the first variation of the action with respect to the 
background solution is zero (for any variation) since 



5C 



6g. 



.fih' 



0; 



6C 



9°,{^°)" 



S^M 



0. 



g°,C^°y 



(2.67) 



Therefore the change in the onshell action is quadratic in the perturbations of the metric 
and matter fields, 



'JS'onsheU = / d'^'^^X 



6C 



5g. 



.fiv 



^9, 



6C 



AgA'P" 



f" + <5$A/ 



A$ 



M 



AgA^'^' 



(2.68) 



For example, in the case of a purely gravitational action, the change in the onshell action 
would be expressed in terms of the Einstein operator acting on the metric perturbation, 
and would be quadratic in the metric perturbation. 

Going back to the case of interest, one wishes to compare the onshell actions for the 
ordinary and superfluid phases, at fixed chemical potential /i and temperature T. The 
metric perturbation is quadratic in the dimensionless ratio fJ,/T, which in turn implies that 
the metric contribution to the onshell action is quartic in n/T. The latter is therefore 
subleading compared to the matter terms: these perturbations are linear in the chemical 
potential, and thus give rise to corrections to the action which are quadratic in fi/T. 

We can illustrate this point by the ordinary phase solution, for which the backreacted 
solution is Reissner-Nordstrom (2.50). As anticipated in the limit of small chemical poten- 
tial the backreaction on the metric is suppressed by factors of /U^2^ ~ n^/T'^. The exact 
expression for the (holographically renormalized) onshell Euclidean^ action in the grand 
canonical ensemble with fj, fixed is 



qE 



/3 



onshell 



2K2.f 



1 + 



7^ 



Vd^i, 



(2.69) 



^The Euclidean continuation of the action is denoted by a superscript. In the static case the Euclidean 
action is obtained by iS — > —S^. 
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with /3 the inverse temperature and Vd-i the regulated volume in the spatial directions. 
The Euclidean onshell action is related to the free energy F as S^ = —f3F. Expanding 
this expression for small ^Zh one obtains 

Csheii = '^r^lV-i (1 + (^ - 2)/^'4 + •••)• (2-70) 

By the logic above we should get the same answer, to this order, by evaluating the action 
for the gauge field in the Schwarzschild background. In other words, the change in the 
action relative to that of the neutral black hole is 

SS = -^ J d^'+'x^^F^^-'F^,, (2-71) 

where g° again denotes the background metric. Evaluating this one indeed obtains 

SSLi^oii = ^^^^(^ - 2)a^', (2.72) 

in agreement with the first correction given above. One could similarly analyse the case 
with fixed charge density p, i.e. the canonical ensemble, by adding the boundary term 

^ldJ:^F,^A'^ (2.73) 

which enforces a boundary condition on n^F^y rather than A^. In this case the onshell 
action is shifted by ^pV^-i, with p being the charge density. 

Let us now turn to condensed phases, beginning with a scalar model in which the matter 
part of the action is 



S=^ I d^'+'xV^g 



-If^'^F^, - |V^ - i<?^^|2 - mVl' 



(2.74) 



Here the complex scalar ^ has mass m? and charge q. For a static solution in which 
At = V{z) and ifj = h{z) with h{z) real, the coupled equations of motion reduce to (2.17). 
The onshell action for the matter fields in the fixed gravitational background g° is 



f^^onshell = Y^ / <f"^^XyJ-g° 



1f^^-F,,-V^{^I:V^^) 



(2.75) 



The scalar term vanishes because 



f dJ^^'ijV f^ip (2.76) 



vanishes both on the horizon (due to ip remaining finite there) and at infinity (since ^ falls 
off fast enough) . Therefore the only term contributing in the onshell action is 



'55'onsheii = Y^ / d^^^xV-5° 






(2.77) 



(together with the appropriate boundary term, if one works instead in the canonical en- 
semble). Comparison of the free energy for the ordinary and condensed phases therefore 
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involves only the kinetic energy of the gauge field. The explicit numerical solutions, see 
appendix D, show that the kinetic energy of the gauge field in the condensed field is such 
that the condensed phase is thermodynaniically favoured. 

For the spin-2 model, substituting (2.15) into the action (2.1) and using the field equations, 
one obtains 









(2.78) 



Working out the boundary term in the second line one finds 

_^/.../^(.^._^), 

However, this boundary term vanishes on the horizon z = Zh since f{zh) = and h{zh) 
is finite (condensed phase) or zero (ordinary phase). The boundary term also vanishes 
at infinity because of the falloff behaviour of the spin-two field: /i ~ z as z — t- while 
f{z) — 7- 1 as z —7- and thus the boundary term scales as a positive power of z as z — t- 0. 
Therefore the onshell action for the spin two model is also given by (2.77), as in the case 
of the scalar condensate. On the one hand this result is not entirely surprising, since the 
equations of motion are the same in the two cases, but the latter did not guarantee that 
the onshell actions would be the same. 

Close to the critical temperature we argued that the condensate scales as e and the back- 
reaction on the metric and gauge fields scales as e^. Since the spin two contribution to 
the action vanishes because of the boundary conditions, the change of the action in the 
condensed phase scales as e^. Recalling that e = (1 — T/TcY'"^ this gives a change in the 
free energy of order (1 — T /TcY, which is what one would expect for a second order phase 
transition. To compute this change in the free energy, the backreaction on the gauge field 
is needed for small chemical potential /u/T ^ 1, while the backreaction on both the gauge 
field and on the metric is needed at finite chemical potential, in which the ordinary phase 
is the charged black hole. 

2.6 Competing phases 

The analysis in the previous sections has an interesting interpretation in terms of competing 
phases. In any top down model, such as those obtained by Kaluza-Klein reductions, there 
will be a multitude of charged fields. Our analysis shows that, for fixed chemical potential, 
scalar and spin two fields having the same mass and charge parameters would condense 
from the ordinary phase at the same temperature. 

If there are no (direct) interactions between the scalar and spin two, then one would expect 
a phase of mixed symmetry, in which both an s-wave condensate and a d-wave condensate 
are present. If one includes direct interactions, and indeed indirect interactions via the 
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metric and gauge field, then the two phases would presumably compete with each other 
and one would need to determine which is indeed the preferred phase. 

In general the phase structure will become very complicated even in a bottom up model 
involving only metric, gauge field, charged scalar and charged spin two field. For example, 
to determine the preferred phase one also needs to take into account that there may be 
qualitatively different solutions of the field equations (as well as the mixed symmetry but 
isotropic phases). 

In a top down model containing towers of charged fields, the phase structure is potentially 
very rich but would also be very difficult to analyse. From (2.59), one can infer that the 
critical temperature behaves as 

(2.80) 



Tc 



r\j 



/i m 

which implies that, for fixed chemical potential, the critical temperature is increased by 
higher charge and lower mass (dimension). This is unsurprising: the lowest dimension op- 
erators become important before operators of higher dimension, and increasing the charge 
corresponds to increasing the coupling to the gauge field. It is interesting to note that, in 
a reduction in which the spin two field is the lightest charged field, it would be the first to 
condense into a superfluid phase. 

3 Kaluza-Klein approach 

In this section we turn to the question of whether and under which conditions the spin two 
model can be obtained from the reduction of a higher dimensional supergravity theory. 

To begin our discussion let us note that massive spin two fields do indeed generically arise 
as Kaluza-Klein modes in reductions leading to anti-de Sitter solutions. To see this let us 
denote the higher dimensional metric as Gmn and consider the most common situation, in 
which the reduction to anti-de Sitter is diagonal over a Sasaki-Einstein manifold. Then the 
higher dimensional metric for the AdS case can be expressed as 

ds^ = Gmndx'^dx'' = g^.dx'^dx-' + g:,dy^dy\ (3.1) 

with Qf^i, the anti-de Sitter metric and g"^^ the Sasaki-Einstein metric, and this metric is 
supported by appropriate fluxes. Such a solution can immediately be generalized to 

ds^ = g^^.dx'^dx" + g^.dy'^dy' (3.2) 

in which g° is the metric on any negative-curvature Einstein manifold M, with the fluxes 
being proportional to the volume forms on the compact and non-compact manifolds. (This 
follows from the fact that Einstein gravity with a negative cosmological constant is always 
a consistent truncation of the higher dimensional theory.) 

The spectrum of the theory is then obtained by considering linearised fluctuations around 
such a background, and then diagonalising the linearised field equations. Since the metric 
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fluctuations are coupled to the flux fluctuations, this diagonalization is rather non trivial 
and will be discussed in detail below for the specific case of compactifications of type IIB to 
five dimensions. It turns out, however, that in the case of the massive spin two fields, the 
diagonalization is in fact immediate: the spin two field arises from the metric fluctuation 

9iiiy = 9% + h^,y. (3.3) 

This fluctuation h^y depends on all coordinates, both x^ and y", and can be expressed 
in terms of the complete basis formed by the harmonics of the Sasaki-Einstein manifold. 
Since it is a scalar from the perspective of the latter manifold, it must be decomposable in 
terms of the scalar harmonics Y^y"') on the compact space as 

K. = Y,hlA^'')Y\y-). (3.4) 

Each such mode results in a spin two field in the lower dimensional space, whose mass 
is given in terms of the eigenvalue of the harmonic on the compact space. (The explicit 
formulae in the case of type IIB reductions will be given below.) The case in which the 
scalar harmonic is trivial corresponds to the massless graviton. Massive spin two fields 
corresponding to Kaluza-Klein gravitons are also charged with respect to the gauge fields 
in the lower dimensional space, since the corresponding scalar harmonics transform in 
non-trivial representations of the isometry group of the compact space. 

Thus the Kaluza-Klein spectrum indeed contains charged massive spin two fields. As 
highlighted earlier, the equations of motion at the linearised level do not however suffice 
to analyse condensate formation. If one works in the neutral black brane background, the 
higher order terms in (2.1), in particular the covariant derivative terms, are necessary to 
generate a condensate solution - without the cubic interaction terms in (2.17) no condensed 
solution exists. However, once one allows for interactions involving the spin two field, then 
one does not expect that the effective action describing the spin two, metric and gauge 
field can involve only a finite number of terms. The question is thus whether from a top 
down perspective one can obtain effective equations of motion equivalent to those obtained 
from (2.1). If one works in the Reissner-Nordstrom background, as in section 2.4, the spin 
two equation is linear, but the gauge field corresponds to a non-diagonal reduction and 
one needs to take into account non-linearities to probe the structure away from the critical 
temperature. 

In this section we show that, at least close to the superfluid transition temperature, the 
bottom up equations can indeed capture the leading terms in the equations of motion. We 
also show that the gyromagnetic ratio g is computable in terms of certain cubic overlaps 
of harmonics on the compact manifold and therefore takes a specific value. Although the 
gyromagnetic ratio does not play any role in the superfluid condensate solution, it is useful 
to derive the value that arises from top down considerations. 
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3.1 Consistent truncations 

A consistent truncation consists of (i) a finite set of lower dimensional fields including the 
metric g^,^ together with other matter fields, denoted schematically by (p^ , which satisfy 
certain equations of motion and (ii) a (non-linear) map relating these lower dimensional 
fields to ten-dimensional supergravity fields {Gmn, Fmnpqr, • • • )• The truncation is said to 
be consistent iff the lower dimensional field equations together with the uplift map give 
solutions of the ten-dimensional equations of motion. 

In the context of sphere, or more generally Sasaki-Einstein, reductions, the consistent 
truncation has traditionally been considered to be a gauged supergravity theory, or a 
consistent subsector of such a theory. More recently novel consistent truncations have 
been found in which certain higher mass fields are retained [37-41]; the corresponding 
operators in the dual CFT are no longer in the stress energy tensor supermultiplet. As 
discussed in [39] the existence of such novel consistent truncations follows from the OPE 
of the dual operators. 

In the case of interest here, however, one does not expect that the lower dimensional action 
can be part of a consistent truncation in which a finite number of fields is retained. Indeed, 
the problems with causality and hyperbolicity of the action (2.1) are indicative that it 
cannot be a consistent truncation. 

On the other hand, the lower dimensional metric and gauge field can by themselves form 
a consistent truncation, with explicit formulas for the reduction being given in the case of 
spherical reductions in [42], see section 5. This fact immediately implies restrictions on the 
couplings between the gauge field and all other Kaluza-Klein fields. Suppose we denote the 
lower dimensional fields by ^'^^ , which consists not just of fields in the lower dimensional 
consistent truncation ip^ , but also all other Kaluza-Klein excitations on the compactification 
manifold. Let us denote the latter by (p , so that ^ = {(j) ,(!)}■ Consistent truncation 
requires that the fields (j) do not source cp i which in terms of the effective lower dimensional 
action implies that there are no couplings involving only one <p 

/"d'^+ix\/^(a"0^)JJap<^^ (3.5) 

where d denotes schematically derivatives, and (n,p) denote integers. 

Conversely any field (p will generically act as a source for many other fields in the reduction. 
Thus for the massive spin two field of interest we would expect that we cannot decouple 
it from the other Kaluza-Klein fields, i.e. a solution of the set of equations following from 
(2.1) would not solve the higher dimensional field equations, as the spin two field would 
source other Kaluza-Klein excitations. 

3.2 Systematic expansion of equations 

Next we review what is known about the interactions between generic Kaluza-Klein fields. 
This question might at first seem intractable since one has an infinite tower of such Kaluza- 
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Klein fields interacting with one another. However, in the literature the interactions have 
been derived using an expansion in the number of fields, for two distinct but related pur- 
poses. 

The first application is in computing correlation functions around a fixed holographic back- 
ground, dual to a given state in the field theory. One then considers all fluctuations to 
be small and in the equations of motion terms involving (^^J^^'^+i are by construction al- 
ways smaller than terms involving ($^^)'^, regardless of the field $*^ under consideration. 
Following the holographic dictionary of [43, 44] quadratic terms give the free action and 
two point functions in the dual theory, with cubic terms needed to compute three point 
functions and so on. In practice only a subset of cubic terms and a small number of quar- 
tic terms have been computed even in maximally super symmetric compactifications, since 
considerable computation is needed. 

A second related use of equations of motion expanded in the number of fields is the follow- 
ing. Suppose one is given an exact solution of the higher dimensional supergravity which is 
asymptotically AdS x SE. Even if this cannot be reduced to a lower-dimensional solution 
of a consistently truncated theory, one can systematically extract holographic data from 
the higher-dimensional solution using the method of Kaluza-Klein holography [45]. The 
basic idea is to express the asymptotic solution near the AdS conformal boundary as a 
perturbation of AdS x SE. Let us focus on a specific example, that of IIB supergravity, 
for which the asymptotic form of the metric and five form should then be 

'^mn — 9mn ~r '^mm ^mnpqr — ^ mnpqr ' Jmnpqri l"^'Oj 

where g'^^ is the AdSxSE metric and F^ is the flux in that background, with hmn being 
the metric fluctuation and fmnpqr being the five form fluctuation. The dilaton ^ = (j)° -\- (j) 
and the axion C = C° + c similarly approach constant values ((/)", C°) near the AdS 
boundary and the asymptotic behaviour of the complex three-form Gmnp = 9mnp is fixed 
similarly. 

Kaluza-Klein holography provides a precise map between the coefficients in the asymptotic 
expansion of the fluctuations (/jmm fmnpqr^ <P, c,gmnp) and the expectation values of chiral 
primary operators in the dual CFT. Operationally this is achieved by expanding these 
fluctuations in the basis of Kaluza-Klein harmonics, and then carrying out a (rather com- 
plicted, non-linear) reduction to obtain the effective five-dimensional action. The crucial 
point is that in obtaining the one point of an operator of dimension A one needs only 
know the non-linear interactions of operators whose dimensions sum to A. For the lowest 
dimension operators, which control the leading behaviour in the IR, this implies that one 
only needs to calculate the lowest order interactions iteratively: cubic, quartic and so on. 

3.3 Spin two equations 

The action (2.1) cannot arise through a consistent truncation but as we discussed above 
that such terms do arise in the expansion of the effective lower dimensional action in terms 
of the number of fields. In writing down the equations of motion (2.17, 2.21), we are using 
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a probe approximation and thus the fields are imphcitly assumed to be small. This would 
suggest that the expansion of the lower dimensional action in terms of the number of fields 
may be valid but this is subtle since the cubic terms in (2.17, 2.21) cannot be viewed as 
negligible relative to the linear terms. If however cubic terms are not subleading to the 
linear terms, then it is not obvious why higher point interactions should be suppressed and 
why other fields in the Kaluza-Klein tower would not be induced. 

Letting $ be a generic Kaluza-Klein field, the issue is that the equations (2.17, 2.21) should 
in general receive additional contributions of the schematic form 

y" + iHj^y _ ^/iV ~ Yl {h^ypq^"^) (3.7) 

n>l,p>l,r>0 

^ ■' ^ ^ ■' ^ ^ n>l,p>l,r>0 

n>l,p>l,r>0 

where £ is a linear differential operator, describing the dynamics of a free non-interacting 
field $. The terms on the right hand side must be at least quadratic in fields and (by the 
argument above) the gauge field cannot act as a source for any of the Kaluza-Klein fields. 
Since the spin two field h is charged, the specific combinations which appear should respect 
charge conservation. Interaction terms will be computed explicitly in the next section; we 
will show that as well as the required cubic terms there are indeed additional terms of the 
above type. 

In the field equation for the generic Kaluza-Klein field $, the r = terms source $, so 
that setting <I> = is no longer a consistent solution. (Terms with r > are less important 
because $ = would always be a consistent solution, although it would not necessarily be 
stable.) In the equations for the spin two field and gauge field, higher point interaction 
terms are not forbidden by any symmetry (apart from charge conversation) and the induced 
Kaluza-Klein fields $ backreact on these equations also. 

Despite the complicated structure of these coupled equations, the terms on the left hand 
side are sufficient to capture the behaviour near the transition temperature. The argument 
can be summarised as follows: following the discussions in appendix C and section 2.4, 
close to the transition temperature the condensed solution admits an expansion 

y = K + eVa + • • • ; /i = e/ii + e^hj, + • • • , (3.9) 

with e = (1 — T/Tc)^ for exponent A. We do not assume that Vq is small compared to the 
temperature; we can take into account the backreaction of the ordinary phase gauge field 
Vo by using the appropriate charged black brane solution; by choosing the U{1) appropri- 
ately we obtain the Reissner-Nordstrom brane, with no other scalar fields induced. In the 
following section we will show that any other Kaluza-Klein field induced by the condensate 
is of order e^ or smaller. This in turn implies that the terms on the right hand side of 
the equation for the spin two are of order e^ or smaller, and thus are subleading to those 
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on the left hand side. In other words, even though other Kaluza-Klein fields are indeed 
generically induced, they do not affect the size of the condensate close to the transition 
temperature. The d-wave phase is not necessarily the most stable phase even close to the 
transition temperature and one would need to analyse the backreaction on other fields to 
explore stability. 

Below the transition temperature, two factors could still suppress contributions to the 
spin two/gauge field equations from other Kaluza-Klein fields. As discussed earlier, large 
charge q can suppress the backreaction of the condensate field. In such a limit, when the 
condensate field is small, but q times the condensate is of order one, only terms involving 
appropriate numbers of powers of q on the righthandside of (3.7) can compete with the 
terms on the left hand side. However, as we will see in the next section, factors of q arise 
from taking derivatives on the compact space. The higher dimensional theory is only two 
derivative and so generically one does not expect higher point interactions to appear with 
more than two powers of q. Making this argument precise, however, would involve the 
complex details of the reduction formulae. Another factor which could suppress the effects 
of the higher point interactions is their finite near horizon behavior. As discussed in the 
appendices, the emergence of a condensed solution at low temperature is rather robust. 
Adding to the spin two equation additional terms which are finite everywhere would be 
unlikely to change the condensed solution qualitatively. 

Analyzing these possibilities in detail, and working out the backreaction of the condensate, 
using the equations (3.7) would be complicated. In section 5 we will discuss the condensed 
phase directly from the higher-dimensional perspective. The backreaction of the condensate 
can be computed directly in the higher dimensional theory and such an approach is likely 
to be easier than working with the dimensionally reduced equations (3.7). For example, the 
Coulomb branch solutions expressed in ten dimensional language are considerably simpler 
than the solutions of the corresponding infinity set of five dimensional equations, analyzed 
in [45]. In that case the complete non-linear solution is best found in ten-dimensional 
language, although note that the reduction to five dimensions is still necessary to extract 
the holographic field theory information, one point functions and so on, see [45]. 

4 Reduction of type IIB supergravity 

In this section we will derive the relevant cubic couplings in the low energy action for 
Sasaki-Einstein compactifications of type IIB supergravity to five dimensions. We will 
determine the gyromagnetic ratio in (2.1) from the top down model and we will show that 
this action indeed captures the leading terms close to the temperature at which the spin 
two condensate forms. 

One might wonder why we choose to work with a five dimensional example, rather than 
one in four dimensions. One reason is for computational convenience as there is more 
literature on couplings of Kaluza-Klein fields in reductions to five dimensions. However, 
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there are conceptual reasons for working with five bulk dimensions rather than four: the 
Coleman-Mermin- Wagner theorem indicates that symmetry breaking in the latter case is a 
relic of the infinite A^ limit. The five dimensional case is stable to finite A^ corrections and 
some d-wave superfluids can in any case be best viewed as having three spatial dimensions, 
for example, certain pnictides. 

An AdS^ X SE^ metric solves the IIB equations supported by a background 5-form flux 
proportional to the volume forms: 

ds^ = g^^dx'^-' + g^^dy'^dy"; (4.1) 

F = 4{r]AdS5 +VSE5), 

where {g^u,9ab) ^^^ *^^ metrics on AdS^ and SE^ respectively, whilst 77 denotes the volume 
form. Our conventions for the type IIB field equations are 

with F = *F. As mentioned in the previous section, this solution can immediately be 
generalized to 

ds^ = g^^^dx'^" + gl.dx'^dx'; (4.3) 

F = 4(r/M5 +VSErJ, 

in which g" is the metric on a negative-curvature Einstein manifold M5. 

4.1 Linearized analysis 

In this section we will discuss the relevant parts of the spectrum of the theory in such a 
background. Let us first consider linearized fluctuations about this background: 

9^ll' ^ 9ij,u "T '^^ll^'■l 9iJ.a ^ /i^aj I4.4j 

9ab = 9ab + ^ab; H^y = h^y + -fl'^j./la 

^mnpq ^ ^mnpq ' ^mnpqi 

where Cmnpq is the 4-form potential and x^ runs over all ten indices {x^, y""). The redefined 
quantity H^^^ takes into account the linearized 5-dimensional Weyl shift needed to obtain 
an Einstein frame metric. 

Here we will not need to consider the fluctuations of other type IIB fields. The latter 
are zero (or in the case of the dilaton and axion, constant) in the background and their 
fiuctuations would only couple to those of the metric and five form via quadratic terms in 
the equations. In other words, the other fields are at least quadratic in the fluctuations 
of interest, {hab,Cmnpq) and the backreaction of these fields into the metric and five form 
equations is then at least cubic order in the fluctuations. 
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Every fluctuation can be expanded in terms of the complete set of harmonics on the compact 
space, for example 

H,. = Y.Hl,{x^)Y\y-), (4.5) 

where Y (y^) are scalar harmonics. Inserting the linearized fluctuations into the field equa- 
tions, and using orthogonality properties of the harmonics, results in linearized equations 
for the five-dimensional fields. In the context of the AdS^ x S^ reduction of type IIB 
the linearized equations were found in [46], and the linear equations given here for AdS^ 
replaced by a negative curvature Einstein manifold are closely analogous. 

There is a gauge dependence, which can be handled by working with gauge invariant 
fluctuations or by fixing a gauge. For computing the spectrum at the linearized level it is 
convenient to impose a de Donder-Maxwell gauge in which 

y^h^ab) = V^haf, = 0; (4.6) 

Va Ty,-. Yy^r- Y/'^n n 

'-■abed — V '-abcfi — ^ '-ab^v — ^ '-a^iup — "J. 

This implies that the complete expansion for the metric and four-form potential fluctuations 
is 

^M- = E^/^-(^)^'(y); V = Y.B';{x)Yl"{y)- 

h(ab) = J2^'Kx)Y;:,){yy, K = Y.^\x)Y\y)- (4.7) 

Cfiupa = / J "/ji^po-l^j^ [y)} Cfj.upa = _/^ "//j/pl^j ^a " (yji 

Cf^uab = Yl ^iS^ (^)^[afe] (y) ; ""f^-bc = ^ 6^" (x)?],,^ V^y/" (y) ; 

Cabcd = Y,^^('')Vab!d^eY^{y). 

Here the harmonics are the scalar harmonics Y^[y), the vector harmonics Y^''[y), the 
traceless symmetric tensor harmonics Y, ', ^ (y) and the antisymmetric tensor harmonics 
yr^y(y). Note that on the 5-sphere these harmonics satisfy 

AY^ = OY^ = -k{k + 4)Y^; (4.8) 

AF/" = (D - 4)yi" = -{k + l){k + 3)y/"; 

^ayl. = ^ayh^^ = ^ayla^^ = Q. (4.9) 

where the symbol A denotes the Hodge-de Rham operator. The indices /, /„, It, la label the 
quantum numbers of scalar, vector, symmetric tensor and antisymmetric tensor harmonics, 
respectively, of given degree k. 
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One can expand the ten-dimensional Ricci tensor Rmn to linear order as 

R^n = R°nn " l^i^x + Oy)hran " ^^mV^/l^ (4.10) 

' rf-'^rn '^rip ~r „ J^^ limp, 

where R^upq is the background value of the Riemann tensor, and Vm is the background 
covariant derivative. O^ and Oy denote the Laplacians on the non-compact space and 
compact space respectively. The background values of the Ricci tensor are R" = —4(7° 
and R°^ = 'ig'^^. The linearized equations can be expressed as 

f(1) = /?(!) A__hklpO pO 771117127713 (4.^^) 

^77171 — -^^77in ' OA mfcmi m2 m,3 n( V^'-^-^/ 

_ __( f Tpomi7n2m'j,7rn , f rpo niim2m3m4\ _ r, 

Qr>\j7nmi7n2m3m4''^7i ~r J 711711111271731714-'^ 771 J — " 

-^mi...m5 = U - /*)mi...m5 + -^h^mi . . .711^ " 5/l[mi-^m2...m5]fc = (4-12) 

where the linearized correction to the Ricci tensor can conveniently be expressed as 

^E = Vfc/lL - ^V™V„,/l|, /iL = \{Vmhi + Vnhi - V''hm7i), (4.13) 

rewriting (4.10). 

Then the linearized field equations can be written down as follows. The linearized Einstein 
equations are 

E^u ■■ ^(n. + ay)H^, - hv,yPH,p + V.VPH^,) + ^V^V^H: - R^.^.h^'^ (4.14) 



\{VnV^Ka + V.V^h^a) - ^<.(n. + n,)/!^ = -^^9% + f 



1 



^2^I7U 



'.,/ . 






i,,i . 



p'^aTfi' V' • 



E^,a ■■ 2 ^°^ ^ °y)hpa - ^{V^VPhap + VaVPHf^p) - —VpVaK + 1^^ f^ aK 

11 8 1 

Eab ■■ T7(°x + ^y)h(^ab) " l^i^ a^Hp + VbVPhap) - —VaVbK + i;'^ a^ hK 






-{VaV'h^bc) + V,V^/i(,,)) + -5°b(nx + ay)h'e - Kcdbh"" 
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Note that Xr^y\ denotes the symmetric traceless part of the tensor X^j^. The hnearized 
five-form equations are 



E, 



flUpCTT 



^d[t.Cupar] = lv^.upar'^"^aCbcde + 2 (^ - ^h^] r]^,p,r] (4.15) 



S, 



pLVpaa • 



(}aCp,upcr ~r '^0\^pCj^pa]a ~ .iVpupcra yprC-bcde + ^C/hC-cdeT) + ^Vpupa'^a.T] 



4! 



E„ 



^d[^C,yp]ab + '2d[aCb]f_ 



1 



Vpupab [p'^cC-deaT ~T ■^C/crCrcde)- 



■'p.upah ■ '~"^[p'^i/p\ab i "^[a^b\p,iip -. ^ ■ 

The hnearized equations for fluctuations of the other type IIB fields decouple from these 
equations. 

Projecting the traceless part of the (//:^) Einstein equation, -Ej-^^,), onto the scalar harmonics 
one obtains 



whilst the trace part of the {^iv) equation, \E'^^ gives 



Y^ = 0, (4.16) 



(4.17) 



16 



' -1 r) ' H' VpCTT ' r (T q 



Y^ =0. 



There are five remaining equations which can be projected onto the scalar harmonics, three 
from the Einstein equations and two from the five- form equations. These are 



E^ 



pa 



IvmL + Ivum' - 4v„7r^ + vy + l^n^'^^hi 



E, 



(ab) ■■ 



K-- 



2 " 15 



pp ' 2 " t^ P 
I 
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V 



"VpCTT 



V(aVf,)y^ = 0; 



Va^' = 0; 

(4.18) 



-□, 



10 150 ^ 5 / 10 ^ '^ 



^y^ 



Y^ =0; 



E^ 



pupar • 
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E, 



l^iupaa 



[bUp.+VpJpayrb']VaY' = 0. 



Y^ =0; 



The equations for {Hf ., H^^,bp^p^,7r^ ,b^) are hence coupled, but one can immediately 
eliminate {H^^,bp^p^) using the second and fifth equations above. Then one obtains cou- 
pled equations for (6^,7r^): 



[(D^ + ay- 32)tt^ + 200^6^] Y^ = 0; 



(4.19) 
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(D^ + ny)b^ - —TT^ 



Y^ =0. 



Let us denote the eigenvalues of the scalar harmonics as OyY^ = —A^Y^; the mass eigen- 
values are then 



m 



Ij = (A^ + 16) - 8^A^ + 4; 



m 



t' 



(A^ + 16) + 8\/A^ + 4, 



(4.20) 
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where the combinations (s , t ) are given by 






(4.21) 



with the masses being defined (n^ — m j)s = (n^ — iTiiji)t = 0. 

The modes {s^,t^) are dual to scalar operators {Ogi,0^i) respectively. These operators 
will acquire no sources or expectation values (at the linearized level) provided that both 
TT^ and b^ are set to zero. Setting these fluctuations to zero however enforces constraints 
on -ff^j,. In particular, from the first two equations in (4.18) one finds that 



{H^y = V^Hi^ = 



(4.22) 



so that the spin two field is both traceless and transverse. These constraints allow (4.16) 
to be simplified, resulting in 



°.-^U-2^(.|p.Hi^«(''-)-AX^)=0. 



(4.23) 



This is indeed precisely the equation (2.5) with iv? = A^. Substituting the Riemann tensor 
in the case of Anti-de Sitter space gives 



with the operator on the left hand side being the linearized Einstein operator. 



(4.24) 



Thus, to summarize, we see that the spectrum necessarily contains transverse traceless spin 
two fields whose mass is given in terms of the eigenvalue of the associated scalar harmonic 
on the compact Sasaki-Einstein space. The presence of such fields does not depend on the 
details of the Sasaki-Einstein space. 

Now let us turn to the vector harmonic sector. The linearized field equations projected 
onto the vector harmonics give 



E, 



a/i ■ 



E, 



(ab) 



E, 






(Am + ^y)Bl^ + \vr^d,bll, + 2Aybjr 
[V^^BJ:] V(,y4^ = 0; 

Kip + r],„'^:V„bi"] Vui;(" = 0. 



i: 



h 



K 



Iv 



0; 



(4.25) 



Ifiup ^ cr^r 



'b] 



where Am denotes the Maxwell operator. The fields bj^^g ^^^ therefore non-dynamical, as 



fiup 



they can be eliminated using the last equation and the remaining equations reduce to 

[(Am + a, - 8)Bl," + AAyb';] Y^ = 0; (4.26) 

[(AM + A,)6^^-4i?^"]y/"=0. 
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These equations can be diagonalized to give two branches: 

Iv — dIv _ M _1_ . /l _1_ A . I h^v. Tn2 



V = ^M - 1 + VI + AZ b';; ml = {4 + AjJ - Ay^lTAT^; (4.27) 



c'; = B'; - (l - v/r+A^j 6^"; m^ = (4 + A,J + 47rTAZ, 

where Aj,y/" = — A/^Y"/" and m? denotes the mass squared of the diagonal combination. 
There is a special case when the vector harmonic corresponds to a Killing vector of the 
compact space, so that A/^ = 8, since then 

a'; =B';-Ah';- m2 = 0, (4.28) 

and ar^ corresponds to gauge fields in five dimensions. Thus we recover the well-known 
fact that gauge fields in the lower dimensional theory are in one-to-one correspondence 
with the Killing vectors of the compact space. 

The field equations associated with the tensor harmonics can also be projected out, result- 
ing in diagonal equations for the scalars ' and the antisymmetric tensors h\2y. From the 
remaining type IIB fields in ten dimensions one obtains a Kaluza-Klein spectrum of scalars, 
massive vector fields and antisymmetric tensors in five dimensions. As argued above, these 
modes are not relevant in what follows and will not be discussed further here; the details 
of this analysis in the AdS case may be found in [46] , 

Let us connect the above discussion to the action (2.1). We have derived diagonal field 
equations for the linearized fluctuations about a background which is the product of a 
negative curvature Einstein manifold and a Sasaki-Einstein manifold. Suppose we switch 
off all such fluctuations, apart from the spin two field H]^^ associated with one specific 
scalar harmonic Y and one gauge field ar^ associated with a Killing vector of the Sasaki- 
Einstein. (The Sasaki-Einstein always admits at least one Killing vector, the Reeb vector.) 
The action for these fluctuations at the quadratic level follows from the reduction of the 
type IIB action over the Sasaki-Einstein manifold and can be computed following the 
analysis in [47] for the case of S^ . 



In particular, the quadratic action for the gauge field a is given by 

SM = -^ f d'xV^^if,,fn, (4.29) 



VSE /■ ,5 /^i 

2kIJ "^''^ ^12^'^^^ 

2 



with VsE the volume of the compact manifold and 1/2k^q the ten-dimensional coupling, 
which are related to the effective coupling in five dimensions given in (2.1) as 



1 _ VsE 



(4.30) 



Here we have assumed that the Killing vectors i^^" are normalized so that 



SE 
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Defining A^ = —=a^ one obtains a canonically normalized gauge field in five dimensions. 
For the corresponding k = \ vector c^ one obtains 

^[«] = -?¥ f d'x./^^{f,Ac)nc) + 48c,cn, (4.32) 

and a canonically normalised massive vector field is obtained by setting C^ = -tsC^. 
The quadratic action for the spin two field H^y = H Y is 

iK^n J 4 2 4 



"10 

(4.33) 

together with the constraints W^ = V ^H^^'^^ = 0. We consider complex scalar harmon- 
ics such that 

/ Y'iY^r = ciVse5'^ , (4.34) 

JSE 

with eigenvalue C\Y = —KY. Then the action can be reduced over the compact manifold 
to give 

S[H] = ^ / d'^^/^c, [-1 V''<V,(i?(^)^^)* - ^mJHl,,{H('^^n* (4.35) 

where rnj = A^. The reason for considering complex scalar harmonics is that they can 
have definite charge under a U{1) symmetry of the compact space. An illustrative example 
would be the following. Consider the five-sphere expressed as 

ds"^ = de^ + sin^ Od^l + cos^ Odip^ . (4.36) 

An example of a complex degree two harmonic satisfying ny(^) = —SY^'^' is 

y(2) = cos^ 06^^^, (4.37) 

which carries definite charge of two under the Killing vector ^^p. The corresponding real 
harmonics 

cos^6lsin(2(^), cos^ 6lcos(2(y9), (4.38) 

do not have definite charge under the symmetry. By rescaling the spin two field in the case 
of a complex harmonic as 

^^.u = -^V^Hf,,, (4.39) 

the total action 

S = S[Af,] + 5[0^,] (4.40) 

gives precisely (2.1) in the free field limit, with the constraints imposed. 

The possible masses of the spin two field and its charge are related to the spectrum of scalar 
harmonics on the compact space. In the case of S^, the scalar harmonics lie in the (0, k, 0) 
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representation of S0{6) and rnj = A^ = k{k + 4), corresponding to a conformal dimension 
A = (4 + /c) with A; > 1 for the spin two field. We will discuss the corresponding range for 
the U{1) charge q below. The corresponding AA = 4 SYM operators are also known 

0^j = Ti (TijX^^^ ■ ■ ■ X^"^^ , (4.41) 

where T^ is the stress energy tensor, X are the 6 scalar fields and the symmetric traceless 
product corresponds to the specific (0, k, 0) spherical harmonic. 

In the case of T^'^ the scalar harmonics are labeled by (j, /, r) of the (SU(2) x SU{2) x f/(l)) 
isometry group and the spin two mass is given by [48-50] 

„2- 



2 - ' •'• ^ \ r 



m'=6l j{j + l) + 1(1 + 1)- j], -=j3 = -/3. (4.42) 

The lowest dimension spin two operators are obtained from harmonics in which I = j = r/2 
for which 

A = 4 + -, (4.43) 

with r a positive integer (and zero for the stress energy tensor). Such fields are in the same 
multiplet as the lowest dimension scalar operators, which have dimension A = 3r/2. The 
effective U{1) charge q of the spin two field however depends on how the gauge field U{1) is 
embedded into the full isometry group; we will discuss q further at the end of this section. 

4.2 Higher order corrections 

Next we turn to the question of the interaction terms in (2.1). Some quadratic contribu- 
tions to the field equations around AdS^ x S^ have been calculated in order to compute 
holographically three-point functions in A/" = 4 SYM, whilst very few cubic contributions 
to the field equations have ever been computed, even around the AdS^ x S^ background. 
Cubic terms involving scalar fields were computed first in [51] and subsequently additional 
cubic terms (and associated three point functions) were obtained in [52]. In the context of 
Kaluza-Klein holography, cubic couplings were discussed in [45, 53, 54]. Here we will build 
on earlier computations to extract the cubic terms of interest in (2.1). 

The corrections to the ten-dimensional field equations up to second order in fluctuations 
can be expressed as 

-t!jmn = J-mni J^mi...m5 = -'mi...m5 (4.44) 

where the quadratic corrections are given by [51, 52] 

Tmi...m5 = -2^lfnii...m5 + ^^[mifm2...mi]k (4.45) 

T^"^) = -i?(2) (-4 4g) 

_i___hkli^s po 7710 mim2m3 _, _i,kisi7^k2S2 po rpo m\m2 

OA I mkmim2m3 ns ' -ta mkik2mim2 nsiS2 

_uks( f rpo ■mim2m-i if rpo mim.2m-i,\ j ^ f f mi.. .7714 

nA \Jmkmim2m,3-'^ns ' J nkmim2m3 -'^ ms ) 1 Q„JmTni...m4Jn 
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where 

R^il = -^kihfhlJ + ^VnikkiVmh'^') + hi^V.hl - hl^hl,, (4.47) 

and h^i was defined in (4.13). 

We now consider the ten-dimensional field equations up to second order in the fluctuations 
of the spin two field and the gauge field, which correspond to the cubic terms in the effective 
action. Let us first note that no term cubic in the spin two field can occur in the action, 
using charge conservation: the spin two field carries a definite U{1) charge and there is 
no cubic combination which is neutral. Similarly the gauge field is Abelian and therefore 
there is no cubic interaction term involving only this field. Therefore the cubic terms of 
interest involve a gauge field, a spin two and its complex conjugate. Looking at (2.1) one 
notes that most of these terms are fixed entirely by gauge invariance, with the parameter 
g being the only free parameter. 

We can extract this cubic coupling from the quadratic corrections to the gauge field equa- 
tion, due to the spin two field -ff(^;/) and its complex conjugate. From (4.25), we see that we 
need to compute the corrections T^^ , TJ^Jpaa and T^Jp^b- However, there are no quadratic 
corrections to the latter two equations from the spin two field. Using the fact that the spin 
two field is transverse and traceless to leading order, we first note that 

42 = -^piH^K^) + \v^{H,pVaH''n - h^Kh^L, (4.48) 



with 



2 



h!^K = \{yahj^)6'^''5K,; (4.49) 






Therefore 



42 = -IVpiHPVaK) + Iv^iH.pVaH^n - \VaH;{V^HP + V^m - V^H, 



1 1 

2 



HP'^VpVaH^, - -VaHP'V^H.p. (4.50) 



(2) (2) (2) 

This correction gives the only contribution to Tap- Tap! = —Rap and thus 

(Am - lQ)Bi^ - 2>2b'^ = -2 tel ,,,, ; (4.51) 



-p """p 

I" _ A R^" 






(Am - 8)b'; - ABJ: = 



where in the first line the correction is projected onto the Killing vector i^a . Note that 
the divergence identities V^i?^" = V^^Tj" = are also uncorrected. Let us define the cubic 
overlap between a Killing vector and the scalar harmonics as 

f K^^^)-Y'Va{Y'r = inCi^ii*VsE, (4.52) 

JSE 
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with Cj^jj* real and n the U{1) charge of the harmonic. One can understand the origin of 
this expression as foUows. Introduce local coordinates for the U{1) such that the Killing 
vector is K = ICd^; the normalization constant /C is then such that the Killing vector is 
normalized as in (4.31). For example, in the case of S^, if one chooses the coordinate (/? in 
(4.36) as the circle direction then the corresponding normalized Killing vector is K = vSc?;^. 
In general, the scalar harmonic can be expressed as Y = e^^'^Y^{9°'), with 9°" the other 
coordinates on the compact manifold, excluding the U{1) circle and hence 

f K^^^)-Y'Va{Y'T = inlCciVsE, (4.53) 

JSE 

where c/ denotes the normalization of the scalar harmonics. The resulting correction to 
the gauge field equation in five dimensions is therefore 



iriK, 



AmA^ = -^m'" {y^.<Ppu - 2Vp<Pf,u) + h.c. (4.54) 



Comparing with (2.3) this gives 



nlC 

71' 



5 = 1. (4.55) 



We can also determine the U{1) charge q and g in five dimensions using the corrections 
to the spin two equation computed in appendix A. This involves computing the gauge 
field strength corrections to the spin two equation at quadratic order, involving one spin 
two field and a vector field. The calculation is somewhat more involved since we need to 
consider the corrections to all seven equations given in (4.16), (4.17) and (4.18) and is 
carried out in appendix A. 

Next we turn to the question of which other fields are sourced by the spin two and gauge 
field. Prom the structure of the equations given above, it is clear that terms quadratic 
in spin two fields also occur in the equations for other Kaluza-Klein fields: the quadratic 
terms can be projected onto different harmonics, and hence different fields. For example, 
from (4.51) one obtains 

(Am - 24)C^ = -^i'l>*r (yf.<Ppu - 2Vp<Pf,,) + h.c. (4.56) 

This implies that the spin two field induces the massive vector field associated with the 
Killing vector. More generally, any cubic interaction which is not explicitly forbidden by 
symmetry is likely to arise. However, the only correction at this order to the gauge field 
equation is that given above - terms involving cp"^ or {(p*)'^ are charged and do not contribute 
to this equation. Close to the transition temperature, when the condensate is of order e, 
the backreaction of such terms on all other Kaluza-Klein fields would be of order e^ . 

The only correction to the spin two equation at this order is discussed in appendix A. 
Terms quadratic in spin two fields can only induce corrections to equations for spin two 
fields with charge 2q, or —2q. Therefore such terms induce spin two fields with other 
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charges, rather than correct the equation for a spin two field of charge \q\. Terms involving 
one gauge field and one spin two field generically source other fields; however, with the 
specific ansatz relevant for the condensed solution, these source terms are shown to vanish 
in appendix A. 

Terms quadratic in gauge fields also generically source other fields: for example, the scalar 
s fields of dimension two are part of the consistent truncation to gauged supergravity 
and they are indeed sourced by the gauge fields. In other words, the backreaction of the 
gauge field even in the ordinary phase sources other fields. By choosing the gauge field 
appropriately, we can restrict to the case in which the backreacted ordinary phase is the 
Reissner-Nordstrom brane, and in this case such scalar fields vanish for the static gauge 
fields of interest. So at finite chemical potential, close to the transition temperature, the 
easiest way to take into account the (finite) effects of the backreaction of the gauge field is 
to work with the uplifted charged black brane which we consider in the next section. One 
may also wish to consider a different choice of ^(1) Killing vector, such that the ordinary 
phase is an R charged black hole, in which scalar fields are turned on. We will discuss this 
possibility later in section 5.3. 

Finally, it is useful to consider the action for charged scalar fields arising in the compact- 
ification. For the fields obtained in (4.20) the effective action has been computed and for 
complex fields of definite charge 

S=^j d'xV^ {-iD,S)iD'^S*) + m2|5p) (4.57) 

^- J d'xV^ {-{D,T){D'^T*) + m'T\T\') , 



2/^2 
where as before the effective charge q appearing in the covariant derivative V^ is 

n/C 
, = -^ (4.58) 

for a field associated with a scalar harmonic of integral charge n. 

4.3 Summary and interpretation 

In this section we have shown that the only terms which contribute to the coupled spin 
two/gauge field equations are those given in (2.16), namely 



(□ - m^)<l)^, = 2R^xpur\ (4.59) 

provided that one (a) imposes 

D^A,<P^''' = = F^pcl)P, (4.60) 

which is implicitly satisfied by the specific static ansatz 

A = V{z)dt; (j)^y = ^h{z) (4.61) 
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and (b) one works with small chemical potential /i/T <C 1 and close to the transition 
temperature. The corrections to the equations which have been dropped are suppressed 
by factors of e^ and (/i/T)^, and the backreaction on the metric and other Kaluza-Klein 
fields is suppressed by the same factors. 

If however the mass and charge parameters are such that condensation does not occur at 
/i/T <^ 1, then it is necessary to drop the restriction to small chemical potential and work 
with the backreacted black brane solution. From the analysis above (and indeed from the 
previous literature), we see that one can indeed not increase q arbitrarily at a given value 
of m? such that the backreaction of the gauge field can be neglected. In other words, 
condensation for spin two fields cannot occur at fj,/T <^ 1. 

More precisely, for the S^, the effective charge q is bounded by the degree of the harmonic 
K as 

\q\ < k, (4.62) 

with equality in the following case. Recalling that the Cartan of S^ is U{1)^, let the scalar 
harmonic have definite integral U{1) charges [rai,n2,n3] with respect to the three U{1) 
factors; group theory implies that |ni| + \n2\ + {n^l < k. Then the harmonic given in (4.37) 
is a degree two harmonic with charges [2, 0, 0] for which q = k = 2 with respect to the first 
L'^(l) factor. Degree k harmonics with charge [k, 0, 0] have q = k with respect to the first 
C/(l) factor. In the context of Reissner-Nordstrom black branes, the relevant harmonics 
are, as we will see in the next section, those of degree k which have equal U{1) charges 
[m, m, m] where m is an integer such that k > 3|?tt.|. Such harmonics give rise to an effective 
charge q = n/v3 with respect to the appropriately normalised diagonal U{1). 

Therefore we find that, for the S^, the charged spin two and scalar fields satisfy 

A^ = {A + k),k>l; As = k,k>2; Ay = 8 + A;, fc > 0, (4.63) 

but in each case \q\ < k. Recall that the critical temperature (2.24) behaves as 

Tc ^ 1 |g| 

fj, TTUc A ' 

with 1 < etc < 2.4. However, for all of the Kaluza-Klein fields 

I < 1, (4.65) 

with the maximum value obtained only for the maximally charged scalar S fields. Therefore 

^c 1 1 , . 

— < < -, (4.66) 

/i TTttc vr 

and therefore the critical temperature can never be reached within the probe approxima- 
tion for the gauge field. In the next section we hence consider spin two fields within the 
rotating D3-brane background corresponding to the uplift of the charged black hole in five 
dimensions. 



(4.64) 
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Even in the backreacted solution, we noted in (2.62) that a condensate only forms at finite 
temperature if the mass to charge ratio is bounded: for D = A, rv? jq^ < 2 at large m. 
This condition gives a strong restriction on which Kaluza-Klein modes can condense since 



-?■ < TV?. 



5 Uplifted solutions 

Close to the transition temperature it is useful to consider the condensate as a linear per- 
turbation around a charged AdS black brane. In this section we will uplift this description 
to a ten-dimensional solution, making use of the nonlinear reduction formulas given in [42]. 
Note that one does not need the complete massless reduction ansatz for type IIB, which is 
still unknown, as it suffices to consider truncations involving only the metric/five form in 
which the gauge fields lie within a Cartan subalgebra of the full gauge group. 

The ten-dimensional metric can be expressed in the form 

ds^ = (g^, + \a^,A^i (ix'^dx" + ds^S'') + -^A^Kadx'^dx'', (5.1) 

where X is a normalised Killing vector in the S^ such that K°^Ka = 1 (using the S^ metric) 
and 

K = d^^+d^^+d^3, (5.2) 

in coordinates for the S^ such that 

3 

ds'iS'') = g^atdx^dx' = Y,{df4 + /^'#?)> (5-3) 

with /ii = sin0, /i2 = cos sin -0 and /X3 = cos 9 cos ijj . 
The five form field strength is F^ + *F^ where 

i^5 = 4e5 + ^ j; d{i^,) A {d(p, + ^) A *5i^, (5.4) 

= 4e5 + l^idK A *5i^) + \Y, di^if) ^A^ *^F, 

i 

where €5 is the volume form on the manifold defined by g^i, and * is the Hodge dual on the 
ten-dimensional space, with *5 being the Hodge dual on the five-dimensional space defined 
by the metric gf^,y. With this reduction ansatz the ten dimensional equations of motion 
are satisfied provided the metric and gauge potential satisfy equations of motion following 
from the lower dimensional action [42] 

^ = 2^ / ^'"^^^ (« + 12 - If^^F,, + -l^e^'^P-'F.^F.^A,^ . (5.5) 

A particular solution of the lower dimensional action, which does not involve the Chern- 
Simons term, is the charged brane solution (2.50), which uplifts to the decoupling limit of 
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the rotating D3-brane solution (with correlated angular momenta). For what follows it is 
useful to note that F^ can in this case be obtained from the potential 

C4 = — jdtAd^x + ^^K AhF, (5.6) 

and thus Fr, only has components Ftzxyw and F^ywab, where the brane spatial coordinates 
are {x,y,w). 

5.1 Linearized perturbations 

In this background solution let us now switch on a generic metric perturbation hmn and 
a five form perturbation fmnpqr- Following (4.10) the perturbation of the ten-dimensional 
Ricci tensor to linear order is 

^mn ~ ~'^'-'il'mn ~ ^^ m^ ni^p ~r ^mpqn^ \^-') 



~r ri (, V J72 V /ipn + V n V /ipm + -ft,^ l^np ~r -tt^ iT'mp) ■ 



where R° denotes the background curvature and □ = V^V"^ is the ten-dimensional oper- 
ator. From (4.11) the linearised Einstein equation is 

^^mn I 9/1 mkmim2in3^ nl \"'"/ 

_( f Tpom\m2m-im,4, . f rpo mim2m-im4\ r, 

r,r>\Jmm,-i_m2m2m4,J^ n i Jnmim2m3m4^rn J — ^• 

yb 

In general, the linearised spectrum will be significantly more complicated than in the case 
where the metric is diagonal. If however we restrict to a metric perturbation hxy{x^,x'^), 
its linearised equation decouples. This follows from the (xy) component of the linearised 
Einstein equation above, which give 

r>(l) _|_ _/^fcZpo pomim2m3 _ PI (^ Q\ 

i.e. hxy decouples from the four form fiuctuations. Moreover, the rest of the linearised 
equations are consistently solved by a transverse hxy, with all other metric fluctuations 
and four form fluctuations vanishing. In other words, just as in the previous section, hxy 
does not couple to and source any other fluctuation. 

Now let us restrict to a static fluctuation which is isotropic along the brane directions, i.e. 
hxy{z,x^). The linearised equation of motion then reduces to 

~ '^'-'hxy + i^xyxy'^ + ^x '^xy + -Z^h ^ yxmi_m2m3^ xy ~ ^' (5. 10) 

where we use R^ = R'y' . This equation needs to be evaluated using the background fields 
given above. Using the background equation of motion this equation can immediately be 
simplified to give 

- ^ahxy + Kyxyh""' = 0. (5.11) 
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We will now show that this equation is equivalent to (2.55). Note that 

^xyxy\9mn) = -^xyxy\9liu) ~ ~~l (5.izj 

and in addition 

det(5'mn) = dei{g^^)dei{gl^) = -j^dei{gl^), (5.13) 

where we have used (2.50). Let the metric perturbation be expressed in terms of a spherical 
scalar harmonic of definite degree k and definite U{\) charge n as 

hxy = \-^I^xy{z)Y\x''), (5.14) 

with 

ny^ = -A;(/c + 4)y^; K-Y^ = inY\ (5.15) 

and c^ denotes the normalization of the harmonic. As discussed at the end of the previous 
section, the three f/(l) charges are equal and integral, i.e. [?tt,, m, m], so that n = 3m with 
|n| < k. 

The equation of motion then reduces to 

((V^ + ig^^)(V^ - iqA'^) - m') <p,y = 2Rly,y<l)-y, (5.16) 

where V^ is the covariant derivative associated with g^y and 

iT? = k{k + A)- 9=^, '^-7^' ^^'^^^ 

This is indeed equivalent to (2.55). The condensate phase is only reached at finite temper- 
ature for small enough masses and large enough charges. 

5.2 Generalizations 

One can immediately generalise the previous analysis to Sasaki-Einstein reductions to five 
dimensions as follows. The metric ansatz is 

ds'^ = g^ydx^dx" + ds\BKE) + (v + A=] , (5-18) 

where r/ denotes the U{1) fibration over the base Kahler-Einstein space Bke and ^ is a 
gauge field on the five-dimensional spacetime with metric g^u- Then the metric 

ds\BKE) + ri^ (5.19) 

is a Sasaki-Einstein metric, provided that the curvature on the base manifold is normalised 
as Rij = 6gij, with rj = dip + P being the Reeb one form on the Sasaki-Einstein and 
Lo = 2^-f being the Kahler form. The base manifold for S^ is CP^ with the Fubini-Study 
metric. The ansatz for the five form is F^ + *F^ where 

Fs = 465 + —=co A hF, (5.20) 

V3 
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where u is the Kahler form on the base space, 65 is the volume form on the five-dimensional 
spacetime defined by the g^y, with *5 being the five-dimensional Hodge dual. The ten- 
dimensional equations of motion are then reproduced by the five-dimensional equations of 
motion following from the action (5.5). 

The reduction from eleven dimensions to four dimensions on a Sasaki-Einstein is very 
similar. The metric ansatz in this case is 

2 



ds^ = g^ydx^'dx^ + ds\BKE) + (v + ^) > (5-21 



where r] denotes the U{1) fibration over the base Kahler-Einstein space Bke (with six 
dimensions) and A is a gauge field on the four-dimensional spacetime with metric 5^,^. 
Then the metric 

ds^{BKE) + V^ (5.22) 

is a Sasaki-Einstein metric, provided that the curvature on the base manifold is normalised 
as Rij = 8gij, with r/ = dip -\- P being the Reeb one form on the Sasaki-Einstein and 
a; = ^dP being the Kahler form. The base manifold for S'^ is CP^ with the Fubini-Study 
metric. The eleven dimensional Einstein equation 

Rran = ^^ (^^M^^n^"^^ " ^/' 9mn) (5.23) 

sets the normalisation of the four form F/^ for which the reduction ansatz can be written 
as 

F4 = 2e4 - a; A *aF, (5.24) 

where uj is the Kahler form on the base space, €4 is the volume form on the four-dimensional 
spacetime defined by the g^^, with =1=4 being the four-dimensional Hodge dual. The eleven- 
dimensional equations of motion are then reproduced by the four-dimensional equations of 
motion following from the action 

5 = ^ y d^xv^ (^R + Q- \f^^ , (5.25) 

provided that one also restricts to solutions which satisfy F f\ F = Q. 

For both eleven and ten dimensional reductions, a specific solution for g^^ and A^ is the 
charged black brane solution given in (2.50). Following the same analysis as in the previous 
section, the following linearised metric perturbations decouple 



+ \^i(t).y{z)Y\x''), (5.26) 



9mn ^ 9mn ■ „ 

where Y {x"") is as before a scalar harmonic on the Sasaki-Einstein. Such metric pertur- 
bations satisfy the equation (5.16), with specific mass and charge parameters, the mass 
for T^'^ being given in (4.42). With the metric parameterisation given above, the charge 
parameter q is proportional to the R charge. 
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5.3 Backreaction 

We now turn to the question of backreaction i.e. given the linearised perturbation can one 
generahse it to find a non-linear solution. The difficulty in finding a non-linear solution 
arises from the fact that the metric perturbation in (5.26) depends not just on the radial 
coordinate of the black brane but also on the coordinates of the Sasaki-Einstein. In general a 
charged scalar harmonic will depend explicitly on all five coordinates of the Sasaki-Einstein. 
For example, in the case of T^'^ one can see this using the metric 

ds^ = ^{d(t) + cos I3idaif + ^ {dl3f + sin^ p^daj) ; (5.27) 

+ -id(l) - COS /32da2) + - idfSo + sin f32da2) ■ 
9 b 

The scalar harmonics can be written in a separable form as [48] 

Yix-) = /,-(/3i)/;(/32)e'^'^"^+'^"^+^^ (5.28) 

where {J3,h) are associated with the SU{2) quantum numbers {j,l) respectively and the 
specific form of the functions fj{/3) is not needed here. Recalling that r/2 = j^, = —I3 we 
note that all harmonics with non-zero R charge depend explicitly on all five coordinates. 
The analysis for the more general case of T^'"^ is very similar. 

To simplify the backreaction it might be interesting to look at the uplift of an R charged 
black branes with scalar fields, rather than Reissner-Nordstrom branes. Using again the 
reduction formulae for S^ given in [42] , let us consider a ten-dimensional metric of the form 

ds^ = A^/^dsl + A"^/2 Ux~^/^de^ + X^^ cos^ e{dip + Af + X^/2 sin^ Odnl) (5.29) 

where 

A = X cos"^ e + X-^/"^ sin^ 9 (5.30) 

and X is a function of the non-compact coordinates on the five-dimensional manifold only. 
When X = 1 and the gauge field A is zero, the metric is a diagonal product on the five- 
dimensional metric with the S^ in coordinates (4.36). With the corresponding ansatz for 
the five form, see [42], the equations of motion are equivalent to those obtained from the 
five dimensional action 



where X = e v2/3x_ This action admits the well-known R charged black brane solution 
ds'^ = -H-'^/^fdt^ + H^'^ {r^dr^ + r^dx ■ dx) ; (5.32) 



A 



with 



sinhp \/x^ 



f = r'H-^; H = l + i^^^^, (5.33) 



42 



with r/i being the horizon location. (To express this metric in a form comparable with that 
of (2.50) we should use a radial coordinate defined as l/z2 = r'^H^'^.) Now x is one of 
the dimension two s scalar fields of the spherical reduction. It is not charged under the 
gauge field retained in the consistent truncation, but since it is associated with a rank two 
harmonic on the sphere it is charged with respect to other U{1) generators in the 50(6) R 
symmetry group and the dimension two operator dual to this field acquires an expectation 
value. 

Now suppose one considers a small metric fluctuation in ten dimensions of the form 

h,y{r,e,^) (5.34) 

which is charged under the Killing vector d^^. From the five dimensional perspective this 
will give rise to a massive charged spin two field. Prom the ten dimensional perspective 
the backreaction is easier than in the previous cases, since the fluctuation depends on only 
two coordinates on the sphere. Moreover condensation of the spin two would necessarily 
compete with that of the scalar field: implicitly there are two order parameters present, 
which may be interesting for modelling the pseudogap region (see conclusions). 

6 Conclusions and outlook 

Charged spin two fields are generic in Sasaki-Einstein reductions and can give rise to con- 
densed d-wave superfiuid phases. Massive spin two fields can never be retained as part 
of a consistent truncation but can be described using a probe approximation, either using 
effective lower dimensional equations or working directly with the ten/eleven dimensional 
supergravity solution. While the condensed solution exists, we were only able to show 
that it was thermodynamically favoured within the probe approximation of the bottom 
up model. Computation of the backreaction of the metric fiuctuation in the higher pic- 
ture would allow us to analyse the thermodynamics and the stability from the top down 
perspective. 

While the focus of this paper was on d-wave phases, our analysis highlights the more 
general fact that Kaluza-Klein modes which are not included in consistent truncations can 
potentially give rise to novel phases in the dual theory. The rule of thumb is that the 
lowest dimension operators control the phase structure but while some such operators are 
retained in consistent truncations to gauged supergravity one needs to take into account 
that other operators of comparable dimension may lie outside these consistent truncations. 

Within our top down approach the masses and charges of spin two fields are not inde- 
pendent and for any given dimension condensed phases only exist at finite temperature 
if the operator charge is large enough. The phenomenology of the condensed phases may 
be explored by analysing the behaviour of probe fermions, for which the masses, charges 
and couplings are also determined by the Kaluza-Klein reduction. The fermion spectral 
functions within the spin two condensed phases will be discussed in our subsequent work 
[55]. 
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As well as the spin two fields discussed here, the top down models of course contain towers 
of other scalar fields, massive vectors etc. The richness of this spectrum could lead to a 
rich landscape of phases of the dual theory, some of which may also not be visible within 
consistent truncations. It would be interesting to explore what phases can be realised and 
to address stability, although the latter is likely to be very difficult given the large number 
of Kaluza-Klein fields. In general one would expect that a number of different phases 
compete with each other. 

As an example of competing phases we highlighted in section 5.3 the case of the spin two 
field competing with a charged scalar field; this would be a relatively simple but potentially 
instructive case to analyse. One of the most mysterious features of high Tc superconductors 
is the pseudogap behaviour in the underdoped regime, just above the critical temperature. 
The physics underlying the pseudogap is not clear, with two of the main theories being 
the existence of preformed Cooper pairs (local rather than long range order) and an exotic 
order parameter competing with superconductivity. The latter scenario is precisely what 
is modelled by including coupled Kaluza-Klein fields and it would be interesting to explore 
this scenario further. 
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A Corrections to spin two equation 

The corrections to the spin two equation at quadratic order, involving one spin two field 
and a vector field, can be used to compute the gyromagnetic coupling. This is somewhat 
involved, as we need to consider the corrections to all seven equations given in (4.16), (4.17) 
and (4.18). These corrections are given by 

Tl^'J = -Rfh (A.l) 

+ \ (V VV"/i.p + V'KaV'Kp) , 
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where we again use the fact that V^hp^ = V^hap = V^hap = /i^ = at the hnearized 
order. The other corrections to the Einstein equations are 

^i? = -Rab = ^^^'V^ (V„ V + Vb V) ; (A.2) 

7^(2) _ _p(2) —h^Pf F° ""^V. 

where the term involving fauarri actually vanishes using the linearized field equations. For 
the five form equations the correction Tpjpar vanishes whilst 

rS.a = 5hlf:p.a]r + Wh^KF^upW, (A.3) 

= ^OhJ^KVpuplrrj- 

In the correction terms we use the leading form for the perturbations, namely 

hpa = BpKa = ^Ka, V = ^M^^' + (H^.Y'T (A.4) 

In explicitly evaluating the quadratic corrections to the spin two equation due to one spin 

(2) 

two field and a gauge field, a considerable simplification occurs: the projection of T^^ onto 
scalar harmonics of the same degree as the spin two field vanishes. For the traceless part 
this follows from the identity 

I y^'"V('^V'')(y^'")*VaK^'") = \{^' - A') inCi^ijVsE. (A.5) 

JSE ^ 

which vanishes when A = A . For the trace part the vanishing of the cubic overlap follows 
from the transverse nature of the vector harmonic. This implies that the corrections to the 
second, third and fourth equations in (4.18) vanish. These uncorrected equations can be 
used to eliminate {T^^,bp^p^) in (4.16) and (4.17). Projecting onto scalar harmonics, the 
correction to (4.16) gives 

^(°. + °y)Ku) + l^i,^u)H^ - V^^V^Hi^^ - Rl\p„\^^H^''' (A.6) 

= ^ [^'(^(^^')P + ^(-^ip - 2Vp^U) - 'iFp^pHi^^ 

From the corrections to (4.17), again projecting onto scalar harmonics one obtains 

1 S inlC 

V^V-HJ,, + -O^H' + -OyH' + 2AH' = —V^A^^HJ,^. (A.7) 

To interpret the terms on the right hand side it is useful to consider the spin two equation 
of motion with a gyromagnetic ratio g, 

(□ - m^)^p,^ = D^<P^ + D^cp^ - Di^^D^)(l) + 2R^p^x^P^ (A.8) 

+gpu [(□ -m"- D)<P - DPcl)p] + igq {Fppr, + F.pO • 
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Expanding this equation up to terms of first order in the gauge field gives 

(□^ - m2)0^^ = V^0^ + V^0^ - V(^,V^)0 + 2R^p^x<t)P^ (A.9) 

+5m. [(°x- -m"- D)(t> - V^p] + igq {F,,p<pP + F.pcpP,) 
-iq{Vp{AP<Pp^) + V^(A^0p,,)) + iqgp^VPA''(l)p^, 



where n^ = V^V^ and in non-linear terms the leading order constraints ^^(j) 



^u 



V^Ap = have been used. 



Comparing the traceless part of this equation with (A. 6) then gives 

1 = ^' 5 = 1, (A.IO) 

in agreement with the values found using the gauge field equation. 
Expanding the constraints given in (2.8) and (2.7) to quadratic order gives 

1(^2 ^ 3)^ _ (^2^ _ v^W^p) = igV^AP^pp; (A.ll) 

^m\m^ + 5)^ = -2iq{g - l)V^F^^V,,/.p^, 

where we use the leading order equation F^ = to simplify the right hand side. Note 
that for g = 1 the second constraint implies (f> = 0, and therefore this constraint remains 
unchanged. Substituting (p = into the other constraint and comparing with (A. 7) again 
gives the same value for q. 

The divergence constraint gives 

m^cPp - Vp^) = -iqm^AP^pp (A. 12) 

-iq{gVpF^,rP + (2 - g)F''PVp(Pp,). 

Manipulating the corrections to the first and last equations in (4.18) gives 

- \/PHp^ + ^V^H = -^^inkFP'^V pH„, 



1 

vIa^" 
1 

73' 



+^=^inkHP''V pF^y (A.13) 



H — -^inkAPRp^. 



which imposing the constraint H = 0, (j) = and recalling that m^ = A^ are consistent 
with the values obtained for q and g. 

Finally let us turn to the corrections to the other equations in this sector, namely the 
equations for the spin one fluctuations a^ and c^; the scalar fields associated the symmetric 
tensor harmonics (f)^^ and the scalar fields associated with scalar harmonics s^ and t^ . 
Projecting the corrections involving one spin two field and one gauge field onto these 
equations we find that these corrections actually vanish, when one uses the explicit ansatz 
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for the spin two field and the gauge field. For example, projecting the first equation of 
(A. 2) onto tensor harmonics gives the correction to the equations for cj)^^. Even if the 
harmonic projection is non-zero, the correction vanishes for the condensate solution since 

H^^'^D^A^ = 0. (A.14) 

The corrections to the other fluctuation equations vanish similarly. 

B Fluctuation equations 

In this appendix we will rewrite the superfluid equations in a form which is useful for 
analysing the critical temperature. Let us start with the following equations of (j){r) and 
ip{r) in terms of a radial coordinate r. 

^"i^) + -^-^'i^) - JJ^'t^ir) = , (B.l) 

These equations are related to the equations in terms of V{C,), h{C,): r = 1/z, (/)(r) = ^(C), 
ip{r) = h{C) and fo{r) := r^(l — 1/r-^), which includes r^ factor compared to f{C) = 1 — C^. 
Note that the last term in (B.l) differs from the previous equations and [56] by a factor of 
2. This is simply related to changing the normalization oiip^ip ^ V^V'- The normalization 
of tp in (B.l) is consistent with that in our main draft. 

The basic idea is to redefine the fields to remove the first derivative terms, which can be 
done by 

(t}{r) = LpQ^j^Lpjr) , %Ij = xo^j^—rT=x{r) , (B.3) 

r 2 T '2 \ Jo 

where ipQ and xo are arbitrary constants. The equations (B.l) and (B.2) then become 

/ + M2(r)vp-A^(r)xV = 0, (B.4) 

X" + M2(r)x + A^(r)(^2^ = 0, (B.5) 



where 



M;- '°-:'f-" . (B.e) 



m 



ip-i)!i , ff r: 



A — ^0 A — '^O (-a 



These equations can be obtained by extremizing the action 



l^'^-lx'^ + lM^^x'-lMy + l^ 



S= I dr{ -ip'^ - -x'2 + -M^x' - -My + -Aip^x' ] , (B.9) 
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where we choose (/?o = XO) so k^p = A^ =: A. 

We can re-express these equations in the C, coordinate using rescaled fields 

y(C) = VoC'^v(c) , h^UoC'^^nO, (B.io) 



so that the equations read 



V" + M^{C)V - kyiOn^V = {) , (B.ll) 

n" + M^{C)n + k'^{c)v'^n = {) (B.12) 



where 



^2 _ p-3)p-l) 



K = - ^2 -^ (B-13) 






2,+^^ .;2,2 h (B-15) 



AV = ^^^ A^ = ^^. (B.16) 



These equations can be obtained by extreniizing the action 

^ = / dC (IV^ - \'H'^ + ^M|/^^ - ^M^V^ + ^AV^?^^ 1 , (B.17) 



where we choose Vq = ^O) so Ay = A-tj =: A. We will use this form of the fluctuation 
equations below. 



C Critical chemical potential (temperature) 



In this appendix we look in detail at solution of the fluctuation equations 



c ef 



V + ^^V - j:-h^V = Q , (C.l) 



^+(7-^ 1^+1 7^-727 1^ = 0, (C.2) 



where recall that 



V = qZHV, h = qh, C=-, ^h = -^, f ={!-(''). (C.3) 

Zh 47rJ 

In the ordinary phase the relevant analytic solution is 

y = /i(i - C^-2) , h = o, (C.4) 
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which is a solution for any /i, with fl defined in terms of the chemical potential fi in (2.32). 
We now analyse the criteria for the existence of other solutions. From the last term of 
(C.2) we may interpret as an effective scalar mass the combination 

mis = m'- jV', (C.5) 

which is nothing but m? + g^^A^. So, if y is big enough, there is a possibility that the effec- 
tive mass of scalar field is sufficiently negative to trigger a BF type instability. Therefore, 
it is natural to start with the form (C.4), plug it into (C.5), and dial up jl until we see the 
instability, i.e. we assume the existence of a perturbative solution to (C.l) and (C.2) near 
the transition temperature of the form [25-27] 

V = Vo + e^V2 + e^i^ + • • • , (C.6) 

h = ehi + e^h3 + e% + --- , (C.7) 

where, Vq = /ic(l — C^'^) ai^d /ic is the critical chemical potential, which gives nonzero hi. 
If hi = 0, then all higher order terms vanish and we end up with the trivial solution (C.4). 

An intuitive argument for the criteria for a condensate solution to exist was given in [4] : the 
effective mass m^g becomes sufficiently negative near the horizon to destabilize the field. 
However, one has to be a little careful about this argument for the following reason. At 
first sight, the effective mass does seem to become increasingly negative near the horizon 
because of the 1// term in (7** or (C.5) which has a pole at the horizon but note that At 
or V has a zero near the horizon and thus there is a vanishing contribution of V to the 
effective mass near the horizon. There is no mass deformation near boundary either due 
to the C^ term in 5**. Indeed, it turns out that a sufficiently negative mass actually arises 
in the intermediate regime, C, ~ 0.5. One exception is the zero temperature case, where g'** 
develops a double pole and g^^V tends to a constant value near horizon. In this case, there 
is a relation between mass and charge for instability obtained by the BF mass violation 
near horizon [4]. 

C.l Method 1: Series expansion and matching 

Working perturbatively close to the transition temperature, we wish to find /ic giving non- 
zero hi: 

( i n _ i\ ■ /f>2 ^2\ 

(C.8) 

Let us solve (C.8) using series expansions near horizon and boundary, which are matched 
at some intermediate point. The relevant series solutions are 
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1.9; 

1.8 ; 



(a) D = 3, 0+ = 2 (b) D = 3, O- = 1 (c) = 4,0+ = 2 

Figure 2. Critical chemical potential, n is the power of terms C" and (1 — (^)" included. The red 
lines are precise values. 

near the boundary and 

hH = ho[l + —il-C) + ^ '-^{l -C? + ---j (CIO) 

near the horizon. 

The matching condition is 

hniU) = hsiU) , h'HiU) = h'BiU) (Cll) 

If ho = 0, then 0_ = 0+ = from (Cll). So we assume /iq 7^ 0. We can choose 
either 0_ = or 0+ = 0. Then we are left with three unknowns 0-)_(0_),/i, Cm and two 
constraints (Cll). i.e. we can have solutions 0+(Cm) and fldCm)- In order to fix the 
ambiguity of Cm, we impose another condition. 

hniCm) = h'siU) (C.12) 

This method is only approximate with a few expanded terms. However, with many terms 
retained, the method becomes very accurate, see figure 2. See figure 3(b) for a comparison 
to other methods. 

C.2 Method 2: Effective scalar in AdS2 or Schrodinger equation 

Next we exploit the heuristic idea of BF mass violation. For this purpose we need to define 
the background for the scalar. First, note that the scalar field <j){z) in AdS^+i reads 

0"_l^<^'_!!j^ = o. (C.13) 

By redefining our field /ii, we can always change the coefficient of h\ of (C.8) to the form 
of the coefficient of (\J in (C.13) with any d. We choose d = 1, which will be interpreted as 
an effective AdS2 scalar equation or Schrodinger equation. 

In terms of the new fields V and 'H, introduced in appendix B 

vb(C) = c'^v(c), h^ = cr^^u(c,), (C14) 
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the leading order and sub-leading order equations read 



P-3)CD-1) 

AC 



v- ^^"7;r^ ^ v=o, (C.15) 



With V, which solves (C.15) with a boundary condition Vo(l) = 0, 

V = k'^{1-Q^-^), (C.17) 

the scalar equation becomes"^ 

'H-Mn = 0, (C.18) 

where 



1^ 



M = - ~^"''^^^'''~\'-^' ^^v-^^-s ; ^c_^g) 



-4(C^TF+ 4g(C-l) +48(-P-^)(^5 + 7£') + ^-^ ^-p^ C~l 



(C.20) 

The BF bound of AdS2 is —1/4. For a large m,^, the effective mass is well above the BF 
bound near the boundary, but saturates the BF bound near the horizon {( — t- 1 — ("), 
regardless of D,m,fi. Therefore, we may expect that the instability can happen by BF 
bound violation near horizon. We define a C dependent mass by 

mls{C) = {l-CfM (C.21) 

where (1 — C)^ corresponds to z^ in (C.13). We find the jlc at a given D and m^ such that 

"^cff (Cc) = -1/4 , m^ff (Cc) = . (C.22) 

The first condition is nothing but the BF bound. The second condition comes from the 
smoothness of m^g((^). These two condition determines flc and <^c uniquely. For example, 
see figure 3(a). At D = 3 and m = 10, the curve with flc = 37.15 and Cc ~ 0.6 satisfy 
(C.22). By repeating for m = 0,2, 5, 10, 15, 20, 25 we can identify corresponding flcS, which 
are plotted as red dots in figure 3(b). 

Let us next consider the limit m? ^ D^. The effective mass then becomes 

^2 / "*- / -2/1 /■D-2\2\ (1 — C) /'ri OQ^ 

^es ~ -pr - /^ (1 - C ) To — ■ (C.23) 



^It can be interpreted as a Schrodinger equation with an inverse square potential. In this case, the 
BF bound violation is interpreted as the existence of negative energy bound state, which is allowed when 
a < —1/4 with a potential — a/r^. 
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q^ ^An 



120 r 




(a) D = 3, 771 = 10, /i = 0, 10, 20, 30, 37.15 (b) D = 3, Line: semi-classical approxima- 
from top tion. Red dots: BF bound analysis. Blue 

dots: numerical values by method 2. 

Figure 3. (a) Effective mass (b) critical chemical potential 



For the purpose of the comparison to a semi-classical approximation, we apply a weaker 
condition 

me'ff(Cc) = 0, m2ff(Cc) = (C.24) 

This condition is equivalent to 



K.ff(Cc)=0, yeff(Cc)=0 



(C.25) 



with 



Vr, 



m 



^/7 



cff 



c 



Mi-C 



D~2^ 



(C.26) 



Note that this equivalence is due to the regularity and positivity of (1 — C)^//^ and the 
reduced condition m'^g = 0. This agrees with the semi-classical analysis in the appendix C 
of [23]. In L> = 3, /ic ~ 3.86771, which is the blue line in Fig. 3(b). The red dots are slightly 
bigger since they are obtained by a stronger conditions mlg{Cc) = —1/4. Interestingly, it 
seems that the semi-classical approximation is quite good for small m, even though the 
condition m,^ ^ D^ is used in the derivation. 
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D Numerical analysis of condensate and comparison to analytic formula 

The purpose of this section is two fold: 1) comparisons between analytic expressions and 
numerical results; 2) checks of the stability and reliability of our Mathematica code. For 
this purpose, we start with the original holographic superconductor model with parameters 
{D = 3,m? = —2, A = 2, K = l,q = 1), for which considerable numerical data has been 
reported, enabling us to cross-check our numerical code with for example [56]. 

D.l Numerical condensate 
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Figure 4. Condensate vs Temperature: Figure (a) agrees with Fig.l in 



D.2 Comparison between numerics and analytic formula 




0.94 0.96 



1.00 T, 



(a) Red curves are by (2.46). (^m ~ (b) Zoom in of (a) near phase transition. 
0.1,0.2,0.26,0.3,0.4,0.5 from top to bot- 
tom. 

Figure 5. Comparison between numerics and analytic formula: the red curves are by (2.46) and 
the blue one is from figure 4(b). 

We note the following: 

• The analytic formula is supposed to be valid only near phase transition, so the devi- 
ations from the numerical result at small temperature are to be expected. However, 
the analytic formula is quite accurate for a much wider regime than would have been 
expected. 
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• The parameter ^^ cannot be fixed in the analytic approximation. It reflects the fact 
that a priori the best matching point is unknown. Therefore, the Cm naay be chosen 
differently depending on the purpose. The parameter was fixed as Cm = 0.26 since it 
gives the best looking fit to the numerics for the widest range. 

• Cm = 0.34 was estimated in [36] not by looking at the plots such as Fig 5, but based 
on other numerical data, T^ ~ 0.118/9^'^, a number borrowed from [56]. However, 
from figure 5, we note that Cm = 0.34 does not yield a good fit. Since there may 
be some numerical uncertainties in the numerical methods between different research 
groups, this discrepancy is arguably not physically important. One possibility is that 
Cm = 0.34 may make the best fit only very near the phase transition. However, if we 
are interested in fermions moving in finite condensate, it is more useful to focus on 
the region a little away from the phase transition point, which motivates the choice 
Cm = 0.26. 

D.3 Comparisons of the bulk fields 

Finally we compare the numerical bulk fields {V, h) to the following analytic expressions 



y{0 



KO 



/^-pC 



D-2 



if < C < Cm, 



vi{l-0 + \{-^ + D- % ) i;i(l - C? if Cm < C < 1- 



h+C' 



if < C < C« 



/^O + ^(1 - C) + ^"(^^-y -"?) (1 - C)^ ifCm<C<l. 



(D.l) 



where /i, p, /i+ are functions of t;i, /iq for a given D, A, C,rn- We choose five points in figure 
6. For each point, there is (t'i,/io) obtained by numerical analysis, which will determine 
(D.l). 




1.0 r,. 



Figure 6. The points where the bulk fields are evaluated. The numbers 1, 2, 3, 4, 5 correspond to 
figure 7. Only the points 1,2,3 are expected to show a good agreement between analytic expressions 
and numerical results. However, we present 4, 5 to show how deviation occurs. (Cm = 0.26 for the 
red.) 
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(a)l 




(b)2 




(c)3 




(d)4 




(e)5 

Figure 7. Comparison between numerical and analytic bulk fields: the red curves are by (D.l) 
and the blue are numerical curves. The analytic curves are not smooth at the matching point. 
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